Cosmological perturbations from stochastic gravity 



Albert Roura 

Theoretical Division, T-8, Los Alamos National Laboratory, M.S. B285, Los Alamos, NM 8754-5 

Enric Verdaguer 

Departament de Fisica Fonamental and Institut de Ciencies del Cosmos, 
Universitat de Barcelona, Av. Diagonal 647, 08028 Barcelona, Spain 

In inflationary cosmological models driven by an inflaton field the origin of the primordial inho- 
mogeneities which are responsible for large scale structure formation are the quantum fluctuations of 
the inflaton field. These are usually computed using the standard theory of cosmological perturba- 
tions, where both the gravitational and the inflaton fields are linearly perturbed and quantized. The 
correlation functions for the primordial metric fluctuations and their power spectrum are then com- 
puted. Here we introduce an alternative procedure for computing the metric correlations based on 
the Einstein-Langevin equation which emerges in the framework of stochastic semiclassical gravity. 
We show that the correlation functions for the metric perturbations that follow from the Einstein- 
Langevin formalism coincide with those obtained with the usual quantization procedures when the 
scalar field perturbations are linearized. This method is explicitly applied to a simple model of 
chaotic inflation consisting of a Robertson- Walker background, which undergoes a quasi-de-Sitter 
expansion, minimally coupled to a free massive quantum scalar field. The technique based on the 
Einstein-Langevin equation can, however, deal naturally with the perturbations of the scalar field 
even beyond the linear approximation, as is actually required in inflationary models which are not 
driven by an inflaton field such as Starobinsky's trace-anomaly driven inflation or when calculating 
corrections due to non-linear quantum effects in the usual inflaton driven models. 

I. INTRODUCTION 

Inflation has become the paradigm for our understanding of the origin of the primordial inhomogeneities which 
are responsible for large scale cosmic structure. The typical inflationary scenario assumes a period of accelerated 
expansion in the early universe, usually driven by a scalar inflaton field, which provides a natural explanation for the 
homogeneity, isotropy and flatness problems of the standard big-bang cosmology [l], 0, H, 0, H[ • The generation of 
structure is explained by the back-reaction effect of the quantum fluctuations of the inflaton field on the gravitational 
field which translate, after quantization, into non-trivial two-point correlation functions of the primordial gravitational 
fluctuations. These correlations give an approximate Harrison-Zeldovich spectrum for large scales frj m. [ToL [Tl| . 
The remarkable success of this scenario to explain the observed anisotropies of the cosmic microwave background 
[3 EH EH EH is today the most compelling reason which supports the inflationary para digm [bollT7|, in spite of some 
interpretational problems such as the transition from quantum to classical fluctuations Il9l. |20L l2lL I22I l23l l24l [25| . 

Semiclassical gravity is a mean field approximation that describes the interaction of quantum matter fields with 
the gravitational field, which is is treated as a classical geometry, and provides a suitable framework for the study of 
macroscopic black holes as well as scenarios in the early universe after the Planck time. In particular it accommodates 
the different inflationary models. The key equation in semiclassical gravity is the semiclassical Einstein equation where 
the expectation value of the stress tensor operator of the quantum matter fields is the source of the spacetime metric. 
In cosmology this is usually assumed to be a spatially homogeneous and isotropic Robertson- Walker spacetime. 
However, since this theory relies only on the expectation value, it completely misses the fluctuations of the stress 
tensor operator. Thus, when the back reaction of the inhomogeneous fluctuations of the inflaton field around the 
homogeneous background are relevant, as in the generation of primordial inhomogeneities, the semiclassical equation 
is insufficient. 

In recent years a stochastic semiclassical gravity, or stochastic gravity, approach has emerged as an extension of 
semiclassical gravity which accounts for the quantum fluctuations of the stress tensor pfil |27| . These fluctuations 
are characterized by the noise kernel, which is defined as the symmetrized two-point quantum correlation function of 
the stress tensor operator. The extension is based on the so-called Einstein-Langevin equation, which is a stochastic 
equation for the linearized gravitational perturbations around a semiclassical background. A Gaussian stochastic 
source with a correlation function determined by the noise kernel is the key ingredient of this equation. From the 
solutions of the Einstein-Langevin equation, the two-point correlation functions for the metric perturbations can be 
obtained. 

Stochastic gravity provides an alternative framework to study the generation of primordial inhomogeneities in 
inflationary models. Besides the interest of the problem in its own right, there are also other reasons that make this 
problem worth discussing from the point of view of stochastic gravity. The Einstein-Langevin equation is not restricted 
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by the use of linearized perturbations of the inflaton field. This may not be very important for inflationary models 
which are driven by an innaton field which takes a non-zero expectation value, because the linear perturbations will 
give the leading contribution; alt houg h the need to consider higher-order corrections (from one-loop contributions) 
has recently been emphasized [28], |29| (see also (3(1). These contributions are in any case important in models such 
as Starobinsky's trace anomaly driven inflation [3l| , which rely on conformally coupled scalar fields with a vanishing 
expectation value. The corresponding Einstein equation is quadratic in these fields and the linear approximation 
becomes trivial. 

In this paper we prove that the usual quantization for linear perturbations of both the metric and the innaton field 
is equivalent to using the Einstein-Langevin equation when the latter is restricted to linearized inflaton perturbations, 
in the sense that the same results for the relevant correlation functions of the metric perturbations are obtained. 

The plan of the paper is the following. In Sec. II a brief description of stochastic gravity is given. This is done in an 
axiomatic way by showing that the semiclassical Einstein equation can be consistently generalized in a perturbative 
way by including a Gaussian stochastic source with vanishing expectation value defined through the noise kernel. The 
dynamical equation for the metric perturbations is the Einstein-Langevin equation. An alternative derivation of this 
equation, reviewed in Appendix C, is based on the influence functional method due to Feynman and Vernon, which 
is g enerally used to describe the dynamics of open quantum systems, i.e., systems interacting with an environment 
[32l [33( 1 . Here the gravitational field plays the role of the system and the quantum matter fields play the role of the 
environment. 

In Sec. Ill we discuss the linearized perturbations around a cosmological Robertson- Walker background coupled 
to a free massive scalar field minimally coupled to the curvature. This corresponds to the simplest model of chaotic 
inflation. The linearized Einstein-Langevin equations is then used to obtain an expression for the correlation function 
of the scalar-type metric perturbations. We concentrate on metric perturbations of scalar type because they are the 
only ones that couple to the inflaton perturbations in the linear approximation. Next, we use the standard linear 
theory of cosmological perturbations, quantize them and derive an expression for the symmetrized quantum correlation 
function of the scalar metric perturbations. This expression is then employed to show the equivalence between this 
correlation function and that derived from the Einstein-Langevin equation. An alternative proof of this equivalence 
is provided in Appendix lEl 

Note that whereas in stochastic gravity the metric is treated as a classical but stochastic field, in the usual ap- 
proach to linear cosmological perturbations both metric and inflaton perturbations are quantized. Nevertheless, the 
correlation functions derived within the Einstein-Langevin approach agree with the symmetrized quantum correla- 
tions, similarly to what happens in simpler open quantum systems [3J|. More specifically, the stochastic correlation 
functions derived from the Einstein-Langevin equation agree with the symmetrized quantum correlation functions of 
the theory of gravity interacting with N matter fields to leading order in 1/N (35[; this was shown in Refs. [36l |37| 
for perturbations around a Minkowski background. In this sense the Einstein-Langevin equation can be regarded as a 
useful intermediary tool to compute the quantum correlations of metric fluctuations in the large N approximation. It 
should also be noted that there are situations for open quantum systems where, for a sufficient degree of environment- 
induced decoherence that guarantees the absence of relevant interference effects, the temporal correlations of some 
actual properties of the system (corresponding to suitably smeared projectors) can be described in terms of classical 
stochastic processes governed by a Langevin equation [38| . In those cases the stochastic correlation functions obtained 
from the Langevin equation also describe such quasi-classical correlations of the system dynamics. 

It is, however, important to stress that when one linearizes with respect to both the scalar metric perturbations 
and the inflaton perturbations, as in the case discussed here, the system cannot be regarded as a true open quantum 
system. The reason is that Fourier modes decouple and the dynamical constraints due to diffeomorphism invariance 
link the metric perturbations of scalar type with the perturbations of the inflaton field so that only one true dynamical 
degree of freedom is left for each Fourier mode. 

In Sec. IV we explicitly compute the correlation functions for the scalar metric perturbations in a simple model 
of chaotic inflation using the Einstein-Langevin approach as described in the previous section. A quasi de Sitter 
expansion for the background is assumed and an almost Harrison- Zeldovich spectrum at large scales is obtained. We 
comment on the different approximations that one is naturally lead to consider within the two approaches. 

Finally in Sec. V we conclude by summarizing our results and briefly discussing the changes that one would 
encounter when using the Einstein-Langevin equation if the inflaton perturbations were treated exactly, that is, 
beyond the linear approximation. We note in particular that the scalar, vectorial and tensorial metric perturbations 
are dynamically related in that case since they all couple to the inflaton perturbations. The metric perturbations can 
then be considered a true open quantum system. 

Throughout the paper we use the (+, +, +) convention of Ref. 39]. We also make use of the abstract index notation 
of Ref. [40(. The Latin indices (a, 6, c...) denote abstract indices, whereas Greek indices are employed whenever a 
particular coordinate system is considered [Latin indices such as (i,j,k...) are used instead when referring only to 
spatial components]. 
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II. STOCHASTIC GRAVITY FORMALISM AND EINSTEIN-L ANGEVIN EQUATION 

A. Stochastic gravity. General formalism 

There are a number of situations, especially in black hole physics and cosmology, in which regarding spacetime 
as classical whereas the remaining matter fields are quantized has proved very fruitful. This is often considered a 
reasonable approximation as long as the typical length scales involved are much larger than the Planck length. A 
first step in that direction is to consider the evolution of quantum matter fields on spacetimes with non-vanishing 
curvature. The consistent formulation of quantum field theory on general globally hyperbolic spacetimes is nowadays 
well established for free fields [4l], [42j , and significant progress has been made for interacting fields as well [HI, 0, [H| . 
Up to this point the quantum matter fields are regarded as test fields evolving on a fixed geometry which is unaffected 
by their presence. 

A second step is to consider the back reaction of quantum matter fields on the spacetime geometry by including 
the expectation value of the stress tensor operator of the quantum fields as a source of the Einstein equation for 
the spacetime geometry, which becomes then the so-called semiclassical Einstein equation. The expectation value of 
the stress tensor is divergent and a non-trivial renormalization procedure is required even for a free field in order to 
preserve general covariance. This can be achieved by using, for instance, dimensional regularization or point splitting 
and introducing suitable local counterterms, which are quadratic in the curvature, in the bare gravitational action. 
Any renormalization method can be used provided Wald's axioms [4l|, [42| are satisfied, since this guarantees equivalent 
results. 

The semiclassical Einstein equation was derived in Ref. [46[ by considering the large N limit of N free scalar fields 
weakly interacting with the gravitational field so that the product of the gravitational coupling constant times the 
number of fields N remains constant as N tends to infinity; see also Ref. [43] for a related result concerning fermions. 
However, the semiclassical Einstein equation is most often introduced in an axiomatic way. The basic aspects of 
this framework, commonly known as semiclassical gravity [42t l48j|. can be summarized as follows. Let us consider 
a manifold M. with a Lorentzian metric g a b which is globally hyperbolic. Let us also consider a linear matter field 
evolving on that manifold. In the Heisenberg picture the scalar field operator ip [g] satisfies the Klein-Gordon equation 

(V a V Q - m 2 ) ${x) = 0, (1) 

where V a means covariant derivative with respect to the metric g a b, and the state of the scalar field, which is 
characterized by a density matrix p [g] , is assumed to be physically acceptable in the sense of Ref. [42j • This means 
that it is of the so-called Hadamard type, so that the expectation value for the stress tensor operator can be consistently 
renormalized. The set (M, g a b, <t> [<?] > P [<?]) constitutes a self-consistent solution of semiclassical gravity if the following 
semiclassical Einstein equation is satisfied: 

G ab [g] + Agab - 2 {aA ab [g] + f3B ab [<?]) = K (f ab [g])' Ien , (2) 

where G a b is the Einstein tensor, (T a b [g])' rcn is the suitably renormalized expectation value of the stress tensor 
operator corresponding to the scalar field operator <f> [g] and a, (3, A and k are renormalized parameters (the prime in 
the expectation value is used to distinguish it from the expectation value introduced below). We considered natural 
units in which h = c = I, and introduced the notation k = 8ttG = 8ir /m~ for the renormalized gravitational coupling 
constant where m p is the Planck mass. The local tensors A a b and B a b are obtained by functionally differentiating 
with respect to the metric terms in the action that correspond to the Lagrangian densities proportional to C abcd C a b c d 
and i? 2 , respectively, where C a b c d and R are the Weyl tensor and the scalar curvature. These terms correspond to 
the finite part of the counterterms introduced in the bare gravitational action to cancel the divergences arising in the 
expectation value of the stress tensor [4l[ . 

From now on, and despite their purely geometric character, we will consider for notational simplicity that the last 
three terms on the left-hand side of Eq. ((2]) have been reabsorbed in the renormalized expectation value of the stress 
tensor operator, which we write now without the prime; this can be done consistently because V a A a b = = V a B a b- 
Taking this into account, the semiclassical Einstein equation becomes 

Gab [g] = K(f a b [ff])ren, (3) 

where we should now keep in mind that the expectation value depends on the renormalized parameters A,a and /3. 

There are, however, situations in which the fluctuations of the stress tensor operator are important |4jJ [5(J [HJ . In 
Refs. [Ill, [53| it was shown that the semiclassical Einstein equation @ could be consistently extended to partially 
account for the fluctuations of the stress tensor operator by introducing a Gaussian stochastic source. More pre- 
cisely, given a self-consistent solution of semiclassical gravity one can introduce the following equation for the metric 
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perturbations h a b around the background metric g a b'- 

Gab [g + h]+A {g ab + hab) - 2 {aA ab [g + k] + l3B ab [g + h}) = K(f ab [g + h])[. cn + K £ab [fl] , (4) 

where the whole equation should be understood to linear order in h a b- Note that throughout this paper indices will 
be raised and lowered using the background metric. The renormalized expectation value is computed with the scalar 
field operator satisfying the Klein-Gordon equation on the perturbed metric g a b + h a b and the Gaussian stochastic 
source £ a & is completely determined by the following correlation functions: 

(Ulg;x))s = o, (5) 

(£ab[g;x)€cd[giy))z = N a bcdix,y) = - ({i ab \g;x) ,tcd\g;y)}) > (6) 

where we used to denote the expectation value with respect to the stochastic classical source £ a b[<?]- The 

operator t a b [g] is defined as t a b [g] = T a b [<?] — (T a b [g]} and the bitensor N a bcd (%> y), which determines the correlation 
function of the stochastic source, is computed using the scalar field operator satisfying the Klein-Gordon equation for 
the background metric g a b- The bitensor N a bcd{x,y) is called the noise kernel, it describes the quantum fluctuations 
of the stress tensor operator and is positive-scmidefinitc. Strictly speaking, the previous definition for the operator 
t a b only makes sense when some kind of regulator is employed since both the operator T a b [g] and the expectation 
value (T a b [<?]) are divergent. However, the operator t a b is finite in the sense that one can compute any matrix element 
of this operator using a regularized version of the two terms that define it, and one finally gets a finite result when 
removing the regulator because the divergences coming from both terms cancel out exactly [53] . Hence, the noise 
kernel requires no renormalization whereas the divergences of the expectation value (T a b [g + h]) appearing in Eq. Q 
are canceled by the counterterms whose finite contribution corresponds to the last three terms on the left-hand side. 
Furthermore, since (£, a b [g])^ = 0, Eq. which is called the Einstein- Langevin equation, reduces to the semiclassical 
Einstein equation for the metric perturbations h a b around the background metric g a b when taking the expectation 
value with respect to the stochastic source ^ab- 1 This framework, in which the metric perturbations are regarded as 
a stochastic process satisfying the Einstein-Langevin equation, is usually referred to as stochastic gravity. 

Similarly to what was done for the semiclassical Einstein equation, we will assume that the last three terms on the 
left-hand side of Eq. (j4]) are reabsorbed in the renormalized expectation value of the stress tensor operator, so that 
the Einstein-Langevin equation will be written from now on as 

G ( ab \g + h] = K(f$ [g + h]) Icn + k U [g] , (7) 

where the superindex (1) means that only terms linear in the metric perturbations h a b are kept. This follows straight- 
forwardly from the fact that Eq. ([4]) was considered only to linear order in h a b (the stochastic source £ a b is regarded 
to be of the same order as h a b) and that the zero order contribution is identically satisfied, since the background 
configuration was assumed to be a solution of semiclassical gravity. 

A necessary condition for the integrability of the Einstein-Langevin equation, via the Bianchi identity, is the 
conservation of the stochastic source. Hence, one must make sure that the stochastic source t; a b [g] is covariantly 
conserved so that Eq. (J7J> is a consistent extension of the semiclassical Einstein equation ©. That the stochastic 
process V a £, a b(x) vanishes is a consequence of the stress tensor conservation on the background metric [52|, [5^ . 
Furthermore, it can also be checked that the Einstein-Langevin equation is compatible with the gauge symmetry 
corresponding to infinitesimal diffeomorphisms. In fact, both the stochastic source and the remaining terms of the 
Einstein-Langevin equation are separately invariant under gauge transformations for the metric perturbations of the 
form h a b — * h a b + V a Ch + Vb^ a corresponding to infinitesimal diffeomorphisms generated by any arbitrary vector field 
C defined on the background spacetime [H, [53| . 

We finish this general introduction to the Einstein-Langevin equation by briefly mentioning that there are derivations 
of the Einstein-Langevin equation in different cosmological settings making use of functional methods [54l . l55l . l56l . l57l . 
[58j or a derivation using arguments based on the renormalization group [5^| . In Appendix [C] we sketch a derivation 
of the Einstein-Langevin equation J71) for the case of a general globally hyperbolic background spacetime using 
the influence functional formalism I5 3| . | 6Q| . The Einstein-Langevin equation has also been applied to the study of 
fluctuations in black hole spacetimes |6ll. l62l. l63j . 



1 Remember that Eq. (J4j should be understood to linear order in h a i,. Therefore, when taking the average over all possible realizations 
of the stochastic source, the equation satisfied by {h a t,)^ coincides with that obtained by linearly perturbing the semiclassical Einstein 
equation J2j- 
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B. Einstein-Langevin equation for cosmological perturbations 

In this paper we will study small perturbations around a Robertson- Walker background when the matter source 
is a minimally-coupled scalar field with a quadratic potential. In fact, this corresponds to the simplest model of 
chaotic inflation with the scalar field playing the role of the inflaton field, but it is sufficient for our purpose of 
illustrating the relationship between the usual treatment of cosmological perturbations and those approaches based 
on the Einstein-Langevin equation within the framework of stochastic gravity. Furthermore, taking into account 
the assumptions made throughout the forthcoming sections, the generalization of our main conclusions and results 
to non-linear potentials should be rather straightforward as long as we keep to quadratic order in the scalar field 
perturbations when considering the potential. 

Recall that the form for the line clement of a general Robertson- Walker metric is 

ds 2 = -dt 2 + a 2 (i) 7ydxW, (8) 

where a (t) is called the scale factor and "ftj is the induced metric for the homogeneous spatial sections, which are 
maximally symmetric hypersurfaces. The line element of the spatial sections can have the three following forms 
r )ijdx l dx^ — {dx 2 + sin 2 xdf2 2 , dr 2 + r 2 d£l 2 , d\ 2 + sinh 2 xcff2 2 } depending on whether the curvature is positive, zero, 
or negative, respectively. In terms of the conformal time coordinate 77 = Jdta _1 (i) the metric (j8]) becomes 

ds 2 = a 2 (77) (-dn 2 + 7^ cfcrW) . (9) 

Before proceeding further it is convenient to introduce the following decomposition for the scalar inflaton field, 
which will be used throughout: 

4,(x) = <t>(r l )+cp(x), (10) 

where 4>(rj), which corresponds to the expectation value (4>[g; x)) of the inflaton field on the background metric, is a 
homogeneous classical-like (as an operator it is proportional to the identity) solution of the Klein-Gordon equation 
which is compatible with the background metric through the semiclassical Einstein equation ([3]). The operator (p(x), 
which will be referred to as the inflaton field perturbations, corresponds to the quantum operator for a minimally- 
coupled massive scalar field whose expectation value vanishes on the background spacetime, i.e., (ip[g;x)} = 0. We 
will consider a Gaussian state for the inflaton field and, thus, for the inflaton field perturbations; see Appendix lAl for 
the definition and the basic properties of pure Gaussian states and the relationship between the state of the inflaton 
field and the inflaton field perturbations. 

It should also be stressed that there are many situations (e.g., in the context of stochastic inflation) in which the 
classical background configuration of the inflaton field will not be homogeneous over the whole spacetime. Nevertheless, 
this will not have observable consequences at present provided that the scale of the inhomogeneities is larger than the 
horizon before the last 60 e-folds of inflation. In fact, when studying models of eternal inflation [64|, [65|, [6y| using the 
formalism of stochastic inflation [67|, the expectation value of the inflaton field is no longer the relevant object. One 
should consider instead the amplitude of a given realization of the inflaton field smeared over scales slightly larger 
than the horizon radius right before the region that had left the self-regenerating regime and would eventually give 
rise to our visible universe underwent the last 60 e-folds of inflation. It has been argued that in those circumstances 
the smeared inflaton field behaves as a classical stochastic process. (This is closely related to the quantum to classical 
transition problem for the inflaton fluctuations 0, [l^, [2(| HH H2, H3, HE]-) If that is the case, one can use a 
particular realization of the smeared inflaton right before the last 60 e-folds of inflation as the classical background 
configuration (t>(i]) and treat it in the same way in which one would have dealt with a quantum expectation value. 

Let us begin by discussing the semiclassical Einstein equation ([3]) for the background metric g a b defined by Eq. ([5]) . 
The right-hand side of Eq. §Ji\) is the properly renormalized expectation value for the stress tensor of the inflaton field 
operator, which satisfies the Klein-Gordon equation |T]) on the background spacetime. 

If we consider the general expression for the stress tensor operator of a minimally-coupled massive scalar field 

fat = Va$V b $ - ^g ab {g cd ^c^d$ + m 2 2 ), (11) 

and use the decomposition of the scalar field introduced in Eq. (fit))) , the expectation value for the stress tensor operator 
can be separated into three different contributions: 

(f ab [g]) Ten = (f ab [g])^ + (f ab [g\)^ + (f ab \g]}™, (12) 

where the subindices <p(j), <p(p and ip<p are used to denote the contributions to the stress tensor operator which are 
respectively quadratic in 4>(r)), linear in both 4>(ji) and (p(x), and quadratic in p(x). The first term depends just 
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on the homogeneous solution 4>[rf), the second term vanishes since it is proportional to ((p[g;x)) and the third term, 
which is completely independent of the homogeneous part <fi(r]), is quadratic in the inflaton field perturbations ip(x) 
and needs renormalization. 

The first term on the right-hand side of Eq. fT2")) will be denoted by T a b = (T a b[g]) see Appendix [Bl for further 
comments on this notation. Taking into account the special form of the Robertson- Walker metric, in the basis 
associated with the conformal time and comoving spatial coordinates these components can be rewritten as 




where primes denote derivatives with respect to the conformal time 77. In this coordinate system the time-time and 
space-space components of a (r))T^ v can be respectively identified with the energy density p(rf) and the isotropic 
pressure p(rj) of a perfect fluid. The components of Eq. (J3]) become then the usual Friedmann equations 

|(V) 2 + mV0 2 ) = 3(W 2 + £ ), (15) 
f (V) 2 - mV0 2 ) = - (2H' + H 2 + e), (16) 

where H = a' /a and e = 0, 1, —1 depending on whether the homogeneous spatial sections of the Robertson- Walker 
geometry are respectively flat, with positive curvature or with negative curvature. 

The third term on the right-hand side of Eq. (fT2")l . (T a b[g])™£, will in turn have a similar structure to that of 
Eqs. (|13[) and (|14[) with diagonal non- vanishing components which can be regarded as corrections Ap(ry) and Ap(j]) to 
the energy density and pressure. This structure is necessary so that the solutions of Eq. are of Robertson- Walker 
type, but there is a family of quantum states of the scalar field which gives rise to such a structure for (Tab ■ 
They can be characterized as follows. Since the Lie derivatives of the six spacelike Killing vectors which characterize 
a Robertson- Walker metric commute with the Klein-Gordon operator satisfying Eq. ([1]) , one can introduce a unitary 
operator which implements at the quantum level the symmetries corresponding to the six Killing vectors and is 
preserved by the dynamical evolution. Consequently, the Hadamard function (the quantum expectation value of the 
anticommutator of the field) employed to compute the renormalized expectation value of the field <p(x) will respect 
the symmetries of the Robertson- Walker geometry provided that one considers a quantum initial state which is kept 
invariant, up to a phase, by the unitary operator associated with those symmetries. Throughout this paper we 
will consider this class of states (spatially homogeneous and isotropic). Nevertheless, being quadratic in the inflaton 
perturbations, which are considered in general to be much smaller during the inflationary period than the homogeneous 
background solution <fi(r]), the contribution from the last term in Eq. (|12[) and, hence, the corrections Ap and Ap, will 
in general be small compared to those from Eqs. (|13[) and (| 1 4j) during the inflationary period. The usual treatments 
which keep to linear order in both the metric perturbations and the inflaton perturbations directly discard them. This 
is actually the situation that we will be interested in here. Therefore, the background solution for the scale factor 
a(rj) is completely determined by Eqs. (|15p and (|16[) without considering the corrections that come from the third 
term on the right-hand side of Eq. (|12p . which is approximated to linear order by (T a b\g]) Ien f=a Tab- 

In addition, either from the conservation of the stress tensor, V a T a b = 0, or by taking the expectation value of 
Eq. (JTJ) , the homogeneous background solution 4>{rf) is seen to satisfy the following Klein-Gordon equation on the 
background Robertson- Walker metric: 

f + 2—<t>' + 7 na 2 <t> = 0. (17) 
a 

Let us now consider the objects which appear in the Einstein-Langevin equation ([7]) and particularize them to the 
case addressed here. The geometric part, i.e., the components of the Einstein tensor for a linear perturbation h a b of 
the metric will be discussed in the next section. The contribution to the expectation value of the stress tensor which 
is linear in the metric perturbation, {T^[g + h]) Ten , can be decomposed according to Eq. (fTU)) as: 

(T^[g + h}) ICD = (fj»\g + h]) H + (fW[g + h])^ + (f£>[g + h])™ (18) 
where the whole inflaton field satisfies now the Klein-Gordon equation on the perturbed metric g a b = 9ab + h a b, 



(V Q V Q -m 2 )0 = O, 



(19) 



7 



and V a means the covariant derivative with respect to g a b- The first term on the right-hand side of Eq. (fT8|) , (T^l 
depends on the scalar field only via the homogeneous background solution (p(rj) , which was already fixed by Eq. (fT7|) 
together with Eqs. dT5j) and ([TO]), and therefore the metric perturbations enter only through the explicit dependence 
of the stress tensor on the metric. Contrary to what happened in Eq. (fTO|) , the second term on the right-hand side 

of Eq. (fT8j) . (X 1 ^)^, no longer vanishes since it is now proportional to (tp[g + h]) and the Klein-Gordon equation 
satisfied by ip[g + h) on the spacetime with the perturbed metric g a b, given by Eq. (119|) . has an inhomogeneous source 
term proportional to the metric perturbation h a b and the homogeneous background solution (f>(rj) which in general 
prevents its expectation value (ip[g + h]) from vanishing. Hence, the only non-vanishing contributions to the second 
term are those which depend implicitly on the metric perturbations through the quantum operator for the inflaton 
perturbations ip[g + h]. 

Finally, the third term, (f$)™, which requires renormalization, will have contributions with either explicit or 
implicit dependence on the metric perturbation. Not only the contributions which depend on the metric perturbations 
explicitly, but also those which depend implicitly via (p[g + h] are ultimately quadratic in the inflaton perturbations 
ip\g] [after solving Eq. (|19|) perturbatively in the metric perturbations]; otherwise they would vanish, as follows from 
the fact that (<p[g]) = 0. 

Similarly to what was said concerning the last term in Eq. (1121) . the last term in Eq. (fT5|) is not taken into account 
by usual approaches to cosmological perturbations, which keep to linear order in the inflaton perturbations as well as 
the metric perturbations. We will not consider these terms either in the next two sections, but some general remarks 
on how to deal with them, and possible implications, will be made in Sec. fVl 

Let us now briefly concentrate on the noise kernel, which accounts for the stress tensor fluctuations and characterizes 
the correlations of the stochastic source £ a (,. It is proportional to ({ia6[ff],icd[ff]}) where t a b = T a b — {T a b), and is 
evaluated on the background metric. Using Eqs. ([TO)) and (fTTj) it can be separated into the following non-vanishing 
terms: 

({iab[g},icd[g}}) = {{Lb{g},i c 4g]})^ + ({t ab [g],t cd [g}}) vv , (20) 

where the first and second terms on the right-hand side are, respectively, quadratic and quartic in the inflaton 
perturbation tp[g]. We used a notation similar to that introduced in Eq. (|12|) since the first term on the right-hand 
side of Eq. (|20[) comes entirely from those contributions to the operators t a b and t c d which are proportional to both 
4>{rj) and <p[g], whereas the last term in Eq. (|2"0)) comes from the contributions to the stress tensor which are quadratic 
in cp[g\. The contribution to the noise kernel which depends on the background homogeneous solution (j>(rj) but not on 
the inflaton perturbation ip vanishes since, being proportional to the identity, the corresponding stress tensor operator 
coincides with its expectation value. The terms linear in ip[g\ also vanish because (<£[g]) = 0. Finally, since we will 
be considering Gaussian quantum states for the inflaton perturbations (see Appendix [S] for a definition and a brief 
description of some basic properties of Gaussian states), the 3-point quantum correlation functions (<p[g]</?[g]<£>[<7]) are 
proportional to the expectation value (<£[<?]) and, therefore, the contributions which are cubic in ip\g] vanish as well. 

It is important to note that both the quadratic and the quartic contributions to the noise kernel are separately 
conserved since both 4>(j]) an d % ) independently satisfy the Klein-Gordon equation (|10[) on the background geom- 
etry; recall that 4>{rf) = {<j)[g\). Due to this fact, we can consistently consider a pair of independent stochastic sources 
£i afc and £ 2q & associated with each term so that £ ah = £i ah + £ 2a f> with (£i a b{x)£i c .d{x')) = ({i ab (x),i cd (x')}), [3], 

(£2 ab(x)& cd(x')) = ({i a b(x),icd(%')}) \g] and (€i a b(z)bcd(%')) = 0. The integrability of the Einstein-Langevin 
equation with any of the two sources is then guaranteed because both sources are separately conserved. 

In the next section we will show that keeping only £i a fc, which can be thought to be of the same order as ip, 
the results obtained using the Einstein-Langevin equation and those from the usual treatments which quantize the 
linearized theory for both the metric and the inflaton perturbations are equivalent. On the other hand, some of the 
main features and consequences of the source £,2ab, which can be regarded as being of quadratic order in <p, will be 
briefly discussed in Sec. fVl and studied in more detail in Ref. [35{ , Of course, when considering the stochastic source 
^2ab, the last term on the right-hand side of Eqs. (TT2"]) and (fTO|) should also be considered since their contribution is 
of the same order as that of £,2ab- 

III. EINSTEIN-LANGEVIN EQUATION FOR LINEARIZED COSMOLOGICAL PERTURBATIONS 

A. Gauge invariant formalism for linearized cosmological perturbations 

Let us consider small metric perturbations around a fixed Robertson- Walker background geometry. It can be shown 
[68| that the most general expression for the components of metric perturbations in some particular coordinate system 
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can be written as 




24> a 2 {r,), 

-(B^ + S^a 2 ^), 

[2 (V>7y - %) + (% + %) + /ly] a 2 (7y), 



(21) 



which depends on ten functions: (p, ij), B and E, the two independent components of each transverse vector Si and 
Fi, and the two independent components of the traceless and transverse symmetric tensor hij. The vectors Si and 
Fi as well as the tensor hij are tangent to the isotropic and homogeneous spatial sections of the Robertson- Walker 
spacetime, but depend in general on the conformal time rj which labels each spatial section. Furthermore, the notation 
\i is used to denote the covariant derivative associated with the metric jij induced on these spatial sections. The 
transversality condition for the vectors and tensor is then written as S 1 ^ = 0, Ft = and h 1 -^ — 0. The global 

factor a 2 (rj) is introduced for later convenience, but could be reabsorbed. The metric perturbations will henceforth 
be treated linearly. 

Four functions describing the metric perturbations are of scalar type, four more are of vector type and finally 
there are two which are of tensor type, according to their transformation properties on the three-dimensional spatial 
sections [68l . l69l . l70| . These types are preserved by time evolution provided that the perturbations of the matter 
sources around the configuration generating the background Robertson- Walker geometry are also treated linearly. 
Those ten functions do not characterize in a unique way non-equivalent perturbed geometries since they may arise 
not only due to real perturbations of the geometry but also to changes of the mapping from the background manifold 
to the perturbed one. Hence, a diffeomorphism generated by a vector field £, considered to be of the same order as 
the metric perturbations, would give an extra contribution C^g a b to the metric perturbation h a b, where g a b is the 
background metric. These local diffeomorphisms do not preserve in general the scalar, vectorial or tensorial nature 
of the metric perturbations. 

There are different approaches to overcome the difficulties derived from this gauge freedom. One approach is to fix 
the gauge [7l| so that further changes on the metric perturbations resulting from coordinate changes, are not allowed. 
This can be achieved by fixing some of the ten functions characterizing the components of the metric perturbations 
either directly specifying some components of h a b or imposing relations between them. A second approach, first 
used by Bardeen [69j], is based on the introduction of so-called gauge-invariant variables which corresponds to using 
linear combinations of those ten functions which remain invariant to linear order under diffeomorphisms generated 
by any vector field £. One can always argue that those gauge-invariant variables coincide with the value taken by 
the functions appearing in Eq. (|2ip (or some linear combination of them) in some particular gauge, as follows from 
the remark that the components of any tensorial object referred to a particular and fixed coordinate system do not 
change when are reexpressed in terms of some new coordinates [72| . 

From now on we will consider spatially flat Robertson- Walker metrics, i.e. jij = Sij in Eq. and concentrate on 
scalar- type metric perturbations. The motivation for the latter is that scalar- type metric perturbations are the only 
ones which couple to matter sources characterized by scalar functions when both metric perturbations and matter 
perturbations (the inflaton perturbations in our case) are treated linearly. Only two true kinematical degrees of 
freedom (i.e., before imposing the Einstein equation) exist for this type of perturbations, in the sense that from 
the four arbitrary functions characterizing scalar metric perturbations, the equivalence classes invariant under local 
diffeomorphism transformations are completely characterized by two arbitrary functions (68l . l73l . W^ . l75j . A particular 
example corresponds to the following two linear combinations of the four functions 4>, tp, B and E, which are invariant 
under local diffeomorphisms: 



These gauge-invariant variables were first introduced by Bardeen [69( with the notation &a = & and = — vp. 

One can also define a gauge-invariant version of the linear perturbations of the Einstein tensor, (C[^)*, which 
depends only on the gauge-invariant functions $ and \P, and is invariant under the same kind of local diffeomorphism 
which preserve the scalar nature of the metric perturbations characterized by $ and '5; see Ref. [zH [zB] for details. 

In fact, the gauge-invariant perturbations (Gj^)„ of the Einstein tensor coincide with the actual components of the 
linear perturbations of the Einstein tensor in the so-called longitudinal gauge, which corresponds to taking E = B = 0. 
Similarly, one could also define a gauge-invariant version of the stress tensor linear perturbations and write a version of 
the Einstein equation for the metric perturbations with both sides explicitly invariant; recall that the whole linearized 
Einstein equation is itself gauge invariant. We will follow an alternative procedure which yields equivalent results. The 



$ = 4>+-[(B 
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idea is to consider the components of the Einstein equation in the longitudinal gauge and notice, as will be explicitly 
shown below, that all the geometric dependence can be written entirely in terms of the gauge invariant variables $ 
and ^> since the only non-vanishing scalar contributions to the metric perturbations in the longitudinal gauge, <f> and 
-0, coincide with $ and 'J. 

In the longitudinal gauge the expression of the perturbed metric for scalar type perturbations on a spatially flat 
Robertson- Walker background in terms of the two gauge-invariant functions $(x) and ^(x) is 

ds 2 = a 2 (ri) [-(1 + 2<f>{x))dif + (1 - 2^{x))5 lJ dx t dx : >] , (22) 

and the components of the linear perturbation of the Einstein tensor are 

G (1) ° = 2a~ 2 (m{H<5> + - V 2 *) , (23) 
G (1) o = 2a- 2 di(H<S> + '$'), (24) 

Sj - a - 2 6 jk d k diD, (25) 



G (1 V = 2a- 2 



(2H' + H 2 ) $ + H& + *" + 2H^' + \^ 2 D 



where 7i — a'(rj)/a(r]), D — <f> — ^, V 2 = S 13 didi and, as we mentioned above, primes denote derivation with respect 
to the conformal time rj. The Einstein equation for the linear perturbations of the metric is 

G w ba=K ( T m_ {T (o )ab)+T{ i )ab y (26) 

where the right-hand side corresponds to those terms which are not included in the background stress tensor 7£ and 
are at most linear in the metric perturbations (some terms from the stress tensor on the background geometry, T^a , 
are present because the scalar field is also perturbed). Note that at the classical level one should simply substitute 
the unperturbed stress tensor T^ for (T^ a b ). However, (T^a) has non-trivial additional components in stochastic 
gravity or when quantizing both the metric perturbations and the scalar field; see the last paragraph of Appendix [B] 
for additional discussion on the notation and related points. Furthermore, from now on we will only consider terms 
which are linear in cither the metric perturbations or the inflaton field perturbations ip, which are both assumed to be 
of the same order. The contribution to the stress tensor of these linear perturbations will be denoted by 5T°. Hence, 
taking all this into account, the components of the Einstein equation for linear scalar perturbations of the metric 
become 

G« a b = KSr h a . (27) 

Let us remember that the expression for the stress tensor of the free massive field ip, which is minimally coupled to 
the spacetime curvature and evolves on the perturbed metric g a b — 9ab + h a b, is 

Tab = V a ^V b (^ - ^g ab (g cd V c pV d <p + mV) . (28) 

The components for the linear perturbations of the stress tensor, ST^, in the basis associated with the conformal time 
and comoving spatial coordinates are then straightforwardly obtained: 

ST ° = a' 2 ($(0') 2 - 0V - m 2 a 2 (f>p) , (29) 
ST* = a- 2 <p'd l ip, (30) 
ST? = a- 2 Si(-$(<P') 2 + (f>'ip' -m 2 a 2 <j)ip) . (31) 

Taking into account the fact that ST/ is diagonal, one can use Eq. (|2T|) [see also Eq. (|25p] with i ^ j to conclude that 
$ and \& are equal [74], [75|, except for a possible homogeneous component (independent of the spatial coordinates), 
which should be included in the background scale factor. Alternatively, the same conclusion can also be reached by 
considering the sum of all the diagonal elements of the ij-components of Eq. I[27p together with Eqs. (|30|) and ([3~T]) 
to substitute tp in terms of $ and which yields V 2 ($ — VP) =0 and hence W = $ provided that they vanish at 
infinity. 

From the Friedmann equations (|15[) and (|16[) for the background solution, we have 

^{4>'Y=H 2 -H\ (32) 
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and we can use this equation to reexpress the terms in Eqs. ([29|) - (|3Tj) which are linear in Substituting these terms 
into the perturbed Einstein equations, given by Eq. (|27[) , we finally get 

V 2 $ - 37i$' - (H' + 2H 2 ) $ = ~ (0V' + mV<^) , (33) 

diiV+H®) = ^4>'dm (34) 

(H' + 2H 2 ) $ + 37i$' + <f>" = | (0V' - m 2 a 2 0</?) ■ (35) 

Similarly, the Klein-Gordon equation for the inflaton perturbations can be obtained by linearizing in both the 
metric perturbations $ and the inflaton perturbations tp the exact Klein-Gordon equation for the whole inflaton field 
4>(rj) + *f{x) on the perturbed geometry, Eq. (|19p . and making use of the fact that the homogeneous background 
solution <j){r}) satisfies the Klein-Gordon equation on the background spacetime, Eq. (fTT)) . It can also be obtained 
from the conservation equation, to linear order in the metric perturbations, of the linearly perturbed stress tensor. 
The result is 

<p" + 2H<p' - VV + m 2 a 2 tp - 4</>'$' + 2m 2 a 2 ^ = 0, (36) 
where we have already taken into account that \& = $. 



B. Einstein-Langevin equation for linearized cosmological perturbations 

Using the results of the previous subsection, the Einstein-Langevin equation can be particularized to the case 
of scalar- type metric perturbations around a Robertson- Walker background geometry with the following result: 

|a 2 ((<5t %+£°) = 3H(H<S> + ^)- V 2 *, (37) 

^ = d^'+H®), (38) 



\ ^ (re 



[2H' + H 2 )<I> + H& + + 2H^' + ^V 2 L» 



4 - \s jk d k diD. (39) 



where we have used Eqs. (|23|) - (f^5|) for the linearized Einstein tensor and we have considered the stress tensor operator 
5T£, which results from keeping terms linear in cither the inflaton perturbations or the metric perturbations. The 
notation (5T^)^, is equivalent to (5T^[g + h]} for the particular case of scalar metric perturbations that we are 
considering. Note, in addition, that all the contributions to (ST^)® are, explicitly or implicitly, proportional to the 
metric perturbations since otherwise they would be proportional to (ip[g]}, which vanishes. In fact, this turns out to be 

important so that {5T^\g + h\) coincides with (T^a[g + h\), which is the object that appears in the Einstein-Langevin 
equation, when one keeps to linear order in the inflaton perturbations, i.e., when only the first two terms on the 
right-hand side of Eq. lfl~8|) are considered. 

The three equations corresponding to the spatial components with equal indices of the Einstein-Langevin equation 
are equivalent due to the symmetries of the Robertson- Walker metric and those of the Gaussian state of the inflaton 
perturbations being considered, which was chosen to be compatible with those symmetries. On the other hand, 
the equation for the spatial components with different indices can be used in a similar way to that of the previous 
subsection in order to show that the gauge invariant functions for the scalar metric perturbations ^(a;) and ^>(x) 
coincide. In this case it is also necessary that the spatial components of the stochastic source £ij(x) with indices i 
vanish identically. Indeed, since £, a b(x) is a Gaussian stochastic process with zero mean, £ij(x) will vanish provided 
that (Ciji x )^cd(y)}^ = 0, which can be argued as follows. The correlation function for £ a f, is defined by the noise 
kernel and, as we are keeping to linear order in the inflaton perturbations, only the first contribution in Eq. (|20p . 
({hj\g],icd\g]}) ^ = ({^iij[g],Si cd [g}}), where 8i ab = 6T ab - (ST a b), should be considered. Finally, ({St i3 [g], Si cd [g]}) 

vanishes for i since 8% j = in that case, as follows from Eq. (|3ip with the inflaton perturbation promoted to a 
quantum operator. Hence, from now on we will take \& = Eq. (|39[) is then trivially satisfied for i and for i =j 
it reduces to (no summation should be understood over the repeated index i) 

| a 2 ((sf;)t, + ej) = {m' + n 2 ) $ + m® 1 + (40) 

It is clear that Eqs. (|3"7|) - (|3"9")) are redundant since we have three equations but only two variables to be determined: the 
function $ characterizing the metric perturbations of scalar type and the expectation value of the quantum operator 



11 



for the inflaton perturbations on the spacetime with the perturbed metric, {(p[g + h]), which will also be denoted in this 
case by (<£) $ . However, despite the apparently excessive number of equations, the system is integrable and solutions 
can be found. This fact is guaranteed by the Bianchi identity provided that the source of the Einstein-Langevin 
equation is conserved. This is indeed the case: the averaged and stochastic sources are separately conserved. On 
the one hand, the conservation of (8T a b)q> is equivalent to the Klein-Gordon equation for the expectation value (<p)$, 
which is completely analogous to Eq. (|36|) : 

(0)1 + 2H (0>; - V 2 + m 2 a 2 (0> 4 - + 2m 2 « 2 # = 0. (41) 

On the other hand, the conservation of the stochastic source is a consequence of the conservation of the noise kernel, 
which in turn relies on the fact that the quantum operator for the inflaton perturbations ip[g] satisfies the Klein-Gordon 
equation on the background spacetime, (V a V a — m 2 ) ip(x) = 0. 

Taking all these considerations into account, the Klein-Gordon equation (|41| can be used to obtain the expectation 
value {tp)$ in terms of We can then easily write the expectation value of the stress tensor linear perturbations 
(8Ta)<s> in terms of $ and use any of the constraint equations, Eq. ([37)) or (|38|) to express $ entirely in terms of 
the stochastic source £ a &; to be specific, in this subsection we will consider Eq. (|3"8f . The spatial derivatives can be 
easily handled by working in Fourier space. Hence, in the rest of this section we will work with Fourier transformed 
expressions in the spatial coordinates. A subindex k will denote the three-dimensional comoving momentum vector k 
that labels each Fourier mode in flat space, i.e., 

*k(v) = I <Pxe- i% -*$>{Ti,x). (42) 

Thus, the Fourier-transformed version of Eq. (|38|) is 

iki(*' k +H* k ) = |a 2 (((<5t < ; ) fc ) $ + (Co)k) , (43) 

where fej is the comoving momentum component associated with the comoving coordinate x l . 

Since (<57^)$ is linear in Eq. (|43p is a first order linear integro-differential equation with an inhomogeneous 
term corresponding to the Oi component of the stochastic source . Therefore, one can always write the solution to 
Eq. gSJ for T) > i] as 

Mv) = $1%) + *$\ti) = + ~ [ V dv'Gi%,v> 2 W) (Co) k iff), (44) 

where is the solution to the homogeneous version of Eq. (|43[) with some given initial conditions at an initial 

(i) . ... (k) 

time ijo, $j/ is a solution of the inhomogeneous equation which vanishes at t]q and G^ et is the retarded propagator 
associated with Eq. (|43|) . The correlation function for the scalar metric perturbation regarded as a solution of the 
stochastic differential equation (|4"3"]l corresponds to 

= (*£\m)*$\m))e + (*£\rh.)*®(m))e + (*$Hm)*$\m))e + (^HmMHm))^ (45) 

where (. . . )^ denotes the average over all possible realizations of the stochastic source, as previously defined. From now 
on we will concentrate solely on the last term, which comes entirely from the solutions of the inhomogeneous equation; 
see Appendix [Dl for a discussion on the role of the initial conditions and the contributions of the homogeneous solution 
to the correlation function. The correlation function has then the following form: 

2 r "n 1 p 

<(® fc i^^M) 

= (|) 2 (G« • (N 0l0z ) kk , ■ (G^f) (Vum), (46) 

where we used the notation A T (ri, if) — A(rf , 77) and A - B = drj A(rj)B(rf), and the factors a 2 (?7i) and a 2 (r]' 2 ) were 
simplified when lowering the spatial indices with the background metric in the last equality. Since we are linearizing in 
the inflaton perturbations, only the first term on the right-hand side of Eq. (f20|) should be considered. The expression 
for the Fourier-transformed version of the noise kernel then becomes 

{N abcd ) kk > M) = \ {{(Si ab ) k (ri), (5i cd ) k , (r/)}) $=0 , (47) 
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where 5i a i, — ST a b — {ST a b), as defined earlier, and (. . .) $=0 is the expectation value for the product of quantum 
operators <5i a b with the field ip evolving on the background metric. We finally obtain the following expression relating 
the correlation function for the metric perturbations and the fluctuations of the stress tensor operator: 

= (f ) 2 \ • ({(si 0i ) k , (^ Oi ) fe ,}> 4=0 • (G^f) ■ (48) 

A detailed example of this kind of computation is given in the next section, where the correlation function of scalar- 
type metric perturbations will be computed for the particular case in which the background solutions for <p(rj) and 
a(rj) correspond to a period of slow-roll inflation. 

We end this subsection by working out the explicit expression for the Oi component of the expectation value (67^) $ . 
From Eq. ((50)) we get 

((«?)*(»?)>• = iha- 2 (ry) {ip k {v)U > (49) 

and everything reduces to compute the expectation value (ipk (??))«• O ne wa Y °f obtaining it is by regarding as 
an external source of the Fourier-transformed version of the linearized Klein-Gordon equation (j4Tj) and solving the 
corresponding inhomogeneous equation perturbatively so that (<pk\g + h]) = (</^[<?]) + (fi^ [g + h]) + 0(h 2 ). The 

expectation value (<f>^ [g]) vanishes, and (ip^[g + h]) is the solution of the inhomogeneous equation with vanishing 
initial conditions, which is proportional to the metric perturbation $fc and can be written as 

= f dv'G^^^'iv'Wkiv') - 2m 2 a 2 (rf) l f>(rf)9 k (rf)j, (50) 

— (k) 

where G^. c l is the Fourier-transformed version of the retarded propagator associated with the Klein-Gordon equation 
I41jl with vanishing initial conditions at 770- One can show that the retarded propagator G^et f° r the Klein-Gordon 



equation (|41[) with Fourier-transformed spatial coordinates is given by 

G&W) = *«V) ([0 k (v),0-k(v')]) o(v - n'\ (51) 
Substituting Eq. (f5"Tj) into Eq. ([5H)l and the result into Eq. (|4U)l one gets the following result for the expectation value 

((ST*)^))* = l k % a- 2 { v )cj ) >(i 1 ) J V drfi ([0 k (r,), 0- k (v')]) «V) tfMW) ~ 2m 2 a 2 ( I) ^(r i > t (r|')) • (52) 

Note that this expression for the expectation value of the stress tensor operator requires no renormalization because 
we linearized with respect to the scalar field perturbations; see Appendix [C] for further comments on this point. 
Furthermore, in Appendix [C] we also show that the expectation value obtained above is in agreement with the general 
expression for the expectation value of the stress tensor which follows from the approach to the Einstein-Langevin 
formalism based on functional methods. 

C. Equivalence with the usual quantization methods 

In this subsection we will show that the result for the correlation function of the metric perturbations obtained in 
the previous subsection using the Einstein-Langevin equation and linearizing in the inflaton perturbations coincides 
with the result which follows from the usual quantization procedures in linear cosmological perturbation theory; see 
for instance Refs. Qlzllzl- 

Let us promote the scalar-type metric perturbations and the inflaton perturbations tp to quantum operators. 
Eqs. (|33[) - (|35[) then become equations for the operators in the Heisenberg picture. In particular, we will concentrate 
on the temporal components of the Einstein equation 

3H$' + 3H 2 $ - V 2 $ = ^a 2 6f Q a \^, (53) 

di(fr + H&) = |a 2 <5T ^, (54) 

where the quantum operator for the inflaton perturbations (p[g + h], on which £7^1,1 depends, satisfies the linearized 
Klein-Gordon equation 

tp" + 2H<p' - V 2 Lp + m 2 a 2 ip - 40'$' + 2m 2 a 2 ^ = 0. (55) 
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The situation is completely analogous to that of the previous subsection except for the fact that the metric perturbation 
is now a genuine quantum operator instead of a stochastic c-number. Thus, taking the Fourier transform for the 
spatial coordinates and proceeding in a similar fashion to the previous subsection, the Klein-Gordon equation (|55|) 
can be solved with the following result: 

MV) = 4%) + / (Wffl&ktf) ~ 2mV(r/)«/,(r/)<M?/)) , (56) 

where ip k ^ (77) is a solution of the Klein-Gordon equation on the background spacetime which contains the entire 

specification of the initial conditions at time r)o and G^t is the retarded propagator with vanishing initial conditions 
associated with the Fourier-transformed version of Eq. (|55[). which coincides with that already obtained for Eq. (|41|) . 

Taking into account Eq. (|56|) . one could use any of the constraint equations (|53|) or ()54|1 to express the quantum 
operator for the metric perturbation $ entirely in terms of the operator for the inflaton perturbations <f>^ [g] (in 
addition to the scalar functions 4>(j]) and a(rj) characterizing the background solution), which satisfies the Klein- 
Gordon equation on the unperturbed geometry. However, before proceeding further it is convenient to discuss some 
useful expressions relating ST a b and its expectation values on both the background spacetime and the perturbed 
geometry. The stress tensor operator 6T a b\$, which is linear in both (p[g + h] and 5>, can actually be written as a 
linear combination of terms proportional to ^°^[<7], the inflaton perturbations on the background metric, and terms 
proportional to $. The latter correspond to terms coming either from the explicit dependence of the stress tensor on 
the metric, which give a local contribution, or from the dependence of ip[g + h] on the metric perturbations according 
to Eq. ([56]) . In fact, since Eq. ([56]) is identical to Eq. (|50|) when substituting the stochastic function $k(v) by the 
operator (?/•), it is clear that all the terms proportional to $ in 8T a b\§ are identical to the terms proportional to $ 
in the operator STab [g + h] considered in the previous subsection. Furthermore, since (tp^lg]) — 0, those terms were 
identical to {5T a b)^, where {ST a b)s, should be understood as the result of replacing $ with <3E> in the expectation value 
(6T a b)$ of the previous subsection. Hence, we have 

8%Lb\& ~ (Sf ab )^ = 6f ab \& =0 = 6fab\&=o ~ (8%ib)&=o, ( 57 ) 

where we used in the last equality the fact that (STab}$= actually vanishes. Eq. (|57|) can be written as 

Sf a b\ $ =Si a b + {5f ab ) $ , (58) 

where 5i a b — ST a b\^ =0 ~ (<%!>) |> = o- ^ should be remarked that taking the expectation value with respect to some 
quantum state of the field ip should be considered with caution here since, due to the constraint equations (|53|) and 
(|54|) . the operators <p and $ are not independent. Thus, strictly speaking, (8T a b)$, should be regarded in this context 
merely as a notation for those terms of 5T~ a b\q, which are proportional to $, in contrast to those proportional to <^( ) [g] . 

Substituting Eq. (|58|) into any of the constraint equations (1531) or (|54|) . one could easily obtain the metric per- 
turbation $ in terms of the operator 6t a b constructed with operators for the inflaton perturbations evolving on the 
background metric. In order to compare in detail with the result of the previous subsection, where the Einstein- 
Langevin equation was used, we will explicitly consider the case in which the constraint equation (|54|1 is used. Having 
substituted 6T a b\§ by St ab + (<5^ofc)$ into the Fourier-transformed version of Eq. (|5"4"|) . we have 

ih^i+Hh) - \ a 2 ((^) fc ) 4 = ^a 2 (6ti) k> (59) 

Taking into account that (ST^)^ is linear in $, one can obtain the following expression for $ in terms of St^ from 
Eq. (55]): 

$ k (ri) = ® i k\ri) + ® ( £ ) (v), (60) 
where $^ h ' ) is a solution of the homogeneous version of Eq. (|59| and 

*%\v) = ~/Va 2 (»/)G&W) W h tf) = I ■ (5t 0i ) k ) (r,), (61) 

where G^l is the retarded Green function and the factor a 2 (77') canceled out when lowering the spatial index i with 
the background metric. In fact, since the homogeneous part of Eq. (|59[) has exactly the same form as that of Eq. (I43p . 
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the retarded Green function G YO l coincides with that of the previous subsection. An analogous remark holds for the 
solution of the homogeneous equation, $| . 

Using Eq. (|6ip and concentrating on the inhomogeneous contribution (a discussion of the homogeneous solution 
and its relationship to the initial conditions is given in Appendix [DJ, the symmetrized two-point quantum correlation 
function for the metric perturbation operator <3E> can be written as 

\({*$Hm), HHm)}) = (f ) 2 \ • ({(&*)„ , (Shi) k ,})^ ■ (G^f) (vum) ■ (62) 

Thus, we can see that the result for the symmetrized quantum correlation function of the metric perturbations coincides 
with the stochastic correlation function (|48|) obtained in the previous subsection using the Einstein-Langevin equation. 

We end this section by making a few remarks concerning the issue of the normalization of cosmological perturbations. 
In principle, one could differentiate Eq. (|54|) with respect to the conformal time (the spatial derivatives can be easily 
eliminated by working in Fourier space) and combine it with Eq. (|53|) to obtain a linear second order differential 
equation for the metric perturbation operator $. However, when trying to quantize a theory beginning with the 
equations of motion instead of an action, one faces a normalization ambiguity which stems from the fact that, 
although any pair of actions that differ in a constant factor yield the same equations of motion (either classical or for 
quantum operators in the Heisenberg picture), their corresponding quantum theories are not completely equivalent. 
In particular, the quantum correlation functions for a given state (e.g. the fundamental state) do not coincide. For 
a linear theory they actually differ by some power of a constant factor which is precisely the square root of the 
proportionality constant between the two actions. This is the reason why in Ref. [74j . when quantizing the theory 
for linear perturbations, the final action was obtained from the original linearized action for a scalar field evolving on 
a metric perturbed around a given background geometry together with the linearized Einstein- Hilbert action for the 
perturbations of that metric. That was done by using the constraint equations to reduce the whole action to that for 
the only true dynamical degrees of freedom. It was precisely in order to avoid the normalization ambiguity explained 
above that such a procedure, which turns out to be rather cumbersome, was used instead of working directly with 
the equations of motion and finding at the end an action which corresponds to the equation of motion for the true 
dynamical degree of freedom. 

On the other hand, the method employed in this section is not affected by such a normalization ambiguity because, 
as can be seen from Eqs. (|54[) and pop , the constraint equation relates the operator for the metric perturbations 
to the operator for the inflaton perturbations on the background metric. The normalization of the latter operator 
is already determined by the usual procedure of quantization starting from the action of a scalar field on a fixed 
spacetime geometry. Hence, the key point was to separate the inflaton perturbation operator tp satisfying the Klein- 
Gordon equation (|55|) into a contribution tfi^ [g] which can be regarded as the inflaton perturbation evolving on the 
fixed background spacetime plus a contribution proportional to the metric perturbation operator <E>, and use then 
the constraint equation to express <& entirely in terms of <^^°- ) [ i g]. In contrast with the approach of Ref. [74| , this 
procedure does not give an explicit expression for the reduced action or even the equation of motion for an isolated 
true dynamical degree of freedom, but it is rather useful (and sufficient) in order to compare with the results obtained 
in the previous subsection by means of the Einstein-Langevin equation. 



IV. PARTICULAR EXAMPLE: COMPUTATION OF THE POWER SPECTRUM FOR LARGE SCALES 

IN A SIMPLE INFLATIONARY MODEL 

In this section we will apply the method developed in Sec. IIIIBI to studying the particular example of metric 
fluctuations induced by the quantum fluctuations of the inflaton field in the context of a simple model of chaotic 
inflation corresponding to a free minimally-coupled massive scalar field. In order to carry out explicit computations, 
we will assume that the Robertson- Walker background geometry is close to the de Sitter geometry. For models with 
exponential inflation, as the one being considered here, this approximation is reasonable during the inflationary period, 
in which the so-called slow-roll parameters controlling the deviations from de Sitter geometry are small, but not for 
later times. One can, nevertheless, obtain useful results from a cosmological point of view because those scales which 
are of cosmological interest at present correspond to scales which left the horizon during the inflationary period. This 
can be understood as follows. On the one hand, the evolution of g ravitational perturbations outside the horizon 
is fairly simple, as can be understood from causality arguments [ill. Wty. and rather independent of the particular 
dynamics of the matter sources. On the other hand, the evolution when the scale reenters the horizon later on during 
the radiation and matter dominated eras has been widely studied using the Newtonian approximation [76|, [77J ■ 

We stress that the results obtained in this section, which are based on the use of the Einstein-Langevin equation, 
are not new. They are basically in agreement with most of the literature based on the simultaneous quantization of 
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gravitational perturbations and inflaton perturbations when both are treated linearly 2 . Of course, this fact ultimately 
follows from the equivalence between both approaches established in Sec. IIII CI (as well as in Appendix [Ef . Thus, the 
purpose of this section is to illustrate with a simple but relevant example how the Einstein-Langevin equation can be 
useful to obtain explicit results concerning cosmological perturbations. 

Let us start by recalling the expression for the Oi component of the Einstein-Langevin equation which was obtained 
in Sec. IIII B I working in Fourier space for the spatial components: 

*fc 4 ($' fe + H<S> k ) = |a 2 (<(«*)*)* + (e ) k ) • (63) 

The expectation value of the linearized stress tensor operator is given by Eq. (|52[) and is non-local in the conformal 
time. In general, this fact makes it difficult to find an analytic expression for the solution of the Einstein-Langevin 
equation. 

One possible approach is to realize that there is a certain linear combination of the different components of the 
Einstein-Langevin equation for which all the contributions from the non-local terms cancel out, as well as those from 
the stochastic source (a detailed proof and discussion is provided in Appendix [E| . In that case the equation that 
one needs to solve, Eq. (|E5[) . is a linear second-order ordinary differential equation. In fact, this equation has the 
same form as one often considered in standard treatments of linearized cosmological perturbations (see Eq. (0.48) in 
Ref. [zl]), and one can take advantage of the existing methods and approximation schemes for solving it. Nevertheless, 
for illustrative purposes we will not follow this approach in this section. We will directly consider Eq. (f6"3")> . neglect its 
non-local part and concentrate on the fluctuating part (neglecting the non-local term is not necessary, but it simplifies 
the problem considerably for a quick calculation). 3 Eq. (|6"3")) then becomes 

2*fc 4 (H<D fe + & k ) ~ ^ (fo) fc , (64) 
m p 

from which we can obtain the metric perturbation $^ in terms of the stochastic source (£oi)fc- We need the retarded 
propagator for the gravitational potential &k, is., the required Green function to solve the inhomogeneous first order 
differential equation (|64p with the appropriate boundary conditions: 

<Fto,rn = -±(eto-rn^+Mrn), (65) 

where f{f},rf) is a homogeneous solution related to the chosen initial conditions. In particular, if we take f{rj,7]') — 
—d(Vo — v') a i r ]')/ a { r l)> we have G]? t (r],ri') = for T) < i] , which gives the stochastic evolution of the metric 
perturbations for r\ > tjq due to the effect of the stochastic source after 770. The correlation function for the metric 
perturbations is then given by the following expression: 

(Mv)^k'(v')h^(^) f d m f d m Gt (v,Vi)G»W,mm 0i (vi)Zk>oi(V2))t. (66) 

The correlation function for the stochastic source is, in turn, connected with the stress tensor fluctuations: 

(tkoMZk'oMi = ~<{*l(»7i),&2)}>^ = (27:) 3 S(k + k')h l k^\ m )<P\ V2 )G^(r ll ,r l2 ) (67) 

where the delta function follows from spatial translational invariance and G*u\t)i, 772) = ({0fc(?7i)> <P— fe(%)}) is the k- 
mode Hadamard function for a free minimally-coupled scalar field which is in a state close to the Bunch-Davies vacuum 



2 More specifically, we will obtain a Harrison-Zeldovich spectrum for the scalar metric perturbations with an amplitude which has the 
right dependence on the parameters of the problem (the Planck mass and the mass of the quadratic inflaton potential). However, our 
simple calculation does not give the right result for the spectral index: it gives a spectral index whose value is exactly one (rather than 
slightly smaller than one). In fact, one can explicitly check that the three main approximations that will be employed in this section 
(namely, neglecting the non-local terms, considering a de Sitter background and computing the quantum correlation function for the 
inflaton field using the massless approximation) all contribute to a comparable deviation from the exact result for the spectral index, 
whereas the correct result is obtained when none of the approximations are made. (All this can be checked by proceeding analogously 
to the calculation in Sec. 8.2.2 of Ref. [7I.) 

3 One can see from our derivations in Sec. IHII that neglecting the non-local term is equivalent to neglecting the terms proportional to the 
metric in Eq. 141t . 



16 



on an almost de Sitter background. The so-called slow-roll parameters account for the fact that the background 
geometry is not exactly that of de Sitter spacetime, for which a(rj) = —1/Hrj with — oo < 77 < 0. 

It is also useful to compute the Hadamard function for a massless field and consider a perturbative expansion in 
terms of the dimensionless parameter m/m p , for which observations seem to imply, as will be seen below, a value of 



the order of 10 6 . Thus, we will consider (771,772) = a{^i)a(rj 2 )G^' (771, 772) = (0|{y fc (?7i), V-k(V2)}\ty such that 
a;t|0) = with yk{">]) = a iji)'Pk{ji) = akUk(i]) + a}_ k u*_ k (ri) and Uk(r]) = (2k)~ 1 / 2 e~ lkv (1 — i/krf) corresponding to the 
positive frequency fc-mode for a massless minimally-coupled scalar field in the Bunch-Davies vacuum state on a de 
Sitter background [4l| . 

The result to lowest order in the mass m of the inflaton field and the slow-roll parameters is: 



(*kfo)<M»/)>< 



m t J m J r,o a{ Jl) a(rf) 

tu-n' [ — ) /,• ''r).:/. : • /,'': 
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(68) 



\{m/m p ); overdots denote here 



where we used the lowest order approximation for <j)(t) during slow-roll: cf>(t) ~ 
derivatives with respect to the physical time t. We considered the effect of the stochastic source after the conformal 
time 770. Notice that the result (|68p is rather independent of the value of 770 provided that it is negative enough, i.e., it 
corresponds to an early enough initial time. This weak dependence on the initial conditions is fairly common in this 
context and can be qualitatively understood as follows: after a sufficient amount of time, the accelerated expansion 
for the quasi-de Sitter spacetime during inflation effectively erases any information about the initial conditions, which 
is redshifted away. The actual result will, therefore, be very close to that for 770 = —00: 



($ fe (77)<M»7')>f - 8? 



771 



fc~ 3 (27r) 3 <5(fc 4- k') cosfc(?7 - 77'). 



(69) 



One remark concerning the massless approximation for the computation of the Hadamard function of the inflaton 
perturbations is needed. It is clear from the equation for the scalar field modes that when one considers scales much 
smaller than the Compton wavelength of the inflaton field, i.e., k/a(rj) 3> m, the effect of the mass term can be 
neglected. On the other hand, for scales larger than the Compton wavelength one could object that the mass term 
should no longer be negligible. However, it can be argued that the mass term can also be neglected for large scales 
provided that the Compton wavelength is much larger than the horizon (the Hubble radius H~ 1 ), i.e., H m. The 
argument goes as follows. For a massless minimally-coupled scalar field in de Sitter spacetime the modes become 
effectively frozen after leaving the horizon: k/a{rj) < H. On the other hand, for a massive scalar field, the modes 
decay approximately like exp(— m 2 At/3H) outside the horizon, but this decay will not be important if m 2 /H 2 is 
small enough. In particular, if 3H 2 /m 2 > 60, the decay factor exp[— (m 2 /3H 2 )H At] will not be too different from 
one for those modes that left the horizon during the last sixty e-folds of inflation (HAt = 60 with At being the 
time between horizon exit and the end of inflation), which includes all the relevant cosmological scales since the 
scale that left the horizon sixty e-folds before the end of inflation corresponds to the size of the visible universe at 
present; any feature with a scale larger than the visible universe appears to us as observationally indistinguishable 
from a homogeneous one. Hence, due to the special behavior of the modes outside the horizon, even when considering 
scales which became much larger than the Compton wavelength before the end of inflation, A:exp(— HAt) < m, it is 
reasonable to approximate a massive scalar field with a massless one as long as H 3> m, which happens to be the case 
in most slow- roll inflationary models and in particular for the simple example considered in this section 

Let us consider the cosmological implications which can be extracted from Eq. (|6"9"|) . especially those related to 
large-scale gravitational fluctuations. These fluctuations are believed to play a crucial role in the generation of the 
large-scale structure and matter distribution observed in our present universe [77| . They are also closely connected 
with the anisotropies in the CMB radiation, which decoupled from matter about 4 x 10 5 years after the Big Bang and 
provides us with very valuable information about the early universe [75|, [7(| . 
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From the analysis of our final result in Eq. (|69|) two main well-known facts can be concluded. First, an almost 
Harrison-Zcldovich scale-invariant spectrum is obtained for large scales. Indeed, for scales clearly outside the horizon 
at the times rj and rj' , i.e., krj, krf <C 1, the right-hand side of Eq. (|69|) becomes proportional to k~ 3 5(k + k') with 
negligible extra dependence on k, ij and rf . Second, since we get ($fc(7y)<i>fc' (?/))£ °c (m/m p ) 2 in agreement with the 
usual results [22|, uMi [zl| , the small value of the CMB anisotropies first detected by COBE imposes a severe bound on 
the gravitational fluctuations, characterized by ($>k(v)^k' (v'))t j which implies the following restriction (fine tuning) 
for the inflaton mass: m/m p ~ 10~ 6 . 

Some comments on the mechanisms considered in earlier related work [78], [7!^ H(J HH which allowed a significant 
relaxation on the fine tuning of that kind of parameters are in order here. In those studies either a self-interacting 
scalar field or a scalar field interacting non-linearly with other fields were considered. The modes of the inflaton 
field corresponding to scales of cosmological interest were regarded as an open quantum system with the environment 
constituted either by the short-wavelength modes in the case of self-interaction or else by other fields interacting with 
the inflaton field. Therefore, one can introduce a stochastic description based on a Langevin equation, as explained in 
Ref. [34[, to study the dynamics of the inflaton field modes. In fact, Langevin-type equations or related stochastic tools 
were employed in the references cited above. Furthermore, it was shown in Ref. [34[ that the validity of the results, 
such as the correlation functions, obtained by those methods is independent of the existence of enough decoherence 
to guarantee the presence of a semiclassical regime for the system dynamics. However, it was also shown that the 
two-point quantum correlation function for the system had two separate contributions (see Eq. (4.9) in Ref. 34]): one 
related to the dispersion of the system's initial state and another one which was proportional to the noise kernel and 
accounted for the fluctuations of the system induced by the interaction with the environment. For natural states in de 
Sitter spacetime such as the Bunch-Davies vacuum (in fact, any reasonable state in de Sitter space tends asymptotically 
to it [82|) the dispersion is proportional to H 2 . This contribution is actually several orders of magnitude lar ger than 
that coming from the term proportional to the noise kernel for the situations considered in Refs. [ill, [7(1, l80l l8ll ] . 
Thus, the two-point quantum correlation function for the inflaton perturbations is dominated by the contribution 
connected to the dispersion of the initial state. This point has been confirmed by a detailed analysis in Ref. [23] ]. 
Moreover, this contribution essentially coincides for the cases of an interacting and a free scalar field. The latter is 
the case being considered throughout this paper and exhibits no noise term for the inflaton dynamics because there is 
no environment for the inflaton perturbations (this should not be confused with the noise kernel for the fluctuations 
of the metric perturbations induced by the quantum fluctuations of the inflaton). One could try to choose the initial 
state, as argued in Ref. (81], so that the contribution from the dispersion of the initial conditions were smaller than 
the fluctuations induced on the modes of the inflaton field, but that would require a great amount of fine-tuning for 
the initial quantum state of each mode, which would become highly unstable due to the large dispersion in momentum 
implied by Heisenberg's uncertainty principle [221 ] and quickly tend to the de Sitter invariant Bunch-Davies vacuum; 
or even have no inflation at all due to the back reaction on the evolution of the background geometry generated by 
such a highly excited state. 



V. DISCUSSION 



In this paper we have studied linearized metric perturbations around a Robertson- Walker background interacting 
with a quantum scalar field and we have shown that, when linearizing the perturbations of the scalar field around 
its background configuration, the Einstein-Langevin equation yields a result for the correlation function of the met- 
ric perturbations equivalent to that obtained in the usual approach based on the linearization and quantization of 
both the metric perturbations and the perturbations of the scalar field around its expectation value. Although for 
the sake of concreteness we have mostly concentrated on the case of a spatially flat Robertson- Walker metric and a 
minimally-coupled scalar field with a quadratic potential, the main result can be generalized rather straightforwardly 
to Robertson- Walker metrics with non-flat homogeneous spatial sections, as well as to general potentials for the scalar 
field and arbitrary coupling to the spacetime curvature. Considering Robertson- Walker metrics with homogeneous 
spatial sections of positive or negative curvature would imply, respectively, the use of three-dimensional spherical or 
hyperbolic harmonics rather than simple Fourier transforms for the spatial coordinates, but that would not substan- 
tially change the procedures and the main conclusions since the basic properties of Fourier transforms employed in the 
text have analogous counterparts for these harmonics [H, [U, [Hj]. On the other hand, the use of a general potential 
should not imply major differences since after all we would linearize with respect to the scalar field perturbations 
around the background configuration. 

In addition, we also provided in Sec. IIVI a particular example illustrating how the Einstein-Langevin equation can 
be used in practice to compute the correlation function at large scales for scalar metric perturbations in cosmological 
inflationary models. In doing so we made use of slightly oversimplified approximations, namely, use of a de Sitter 
background geometry, calculation of the Hadamard function for a massless scalar field and neglecting a non-local 
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term, because the result has already been computed in a number of references, see for instance Refs. [!, [7^, l75j . 
and our primary concern was simply to show how the Einstein-Langevin equation can be used to obtain an explicit 
result for cosmological perturbations. 

Throughout the article we have concentrated on scalar-type metric perturbations. The reason for this is that, when 
linearizing with respect to both the metric perturbations and the scalar field perturbations, the vectorial and tensorial 
metric perturbations decouple from the matter scalar field. In this case the metric perturbations do not constitute 
a true open system since the dynamics of the scalar and vectorial perturbations is completely constrained by the 
temporal components of the Einstein equation; in fact, the vectorial ones actually turn out to vanish and the scalar 
ones cannot be regarded as a degree of freedom independent of the scalar field perturbations. Moreover, the only 
true dynamical degrees of freedom, the two tensorial ones, do not couple to the matter field. On the other hand, 
an even more interesting situation corresponds to the case in which the scalar field is treated exactly, at least for 
quadratic potentials. Then the scalar field also couples to the metric perturbations of tensorial type and the metric 
perturbations become a true open system with the scalar field corresponding to the environment. 

The main features that would characterize an exact treatment of the scalar field perturbations interacting with 
the metric perturbations around a Robertson- Walker background as compared to the case addressed in this paper 
are the following. First, the three types of metric perturbations couple to the perturbations of the scalar field, as 
already mentioned above. Second, the corresponding Einstein-Langevin equation for the linear metric perturbations 
will explicitly couple the scalar and tensorial metric perturbations. Furthermore, although the Fourier modes (with 
respect to the spatial coordinates) for the metric perturbations will still decouple in the Einstein-Langevin equation, 
any given mode of the noise and dissipation kernels will get contributions from an infinite number of Fourier modes 
of the scalar field perturbations (see Ref. [86| for an explicit calculation of the noise kernel for a massless and 
minimally-coupled scalar field in de Sitter). This fact will imply, in addition, the need to properly renormalize the 
ultraviolet divergences arising in the dissipation kernel, which actually correspond to the divergences associated with 
the expectation value of the stress tensor operator of the quantum matter field evolving on the perturbed geometry. 

The importance of considering corrections due to one-loop contributions from scalar field perturbations, beyond 
the tree level of the linear cosmological perturbation theory, has recently been emphasized [H, [2t| . In the present 
context this means treating the scalar field perturbations exactly in the Einstein-Langevin equation. Furthermore, 
in Ref. [34j it was explained how a stochastic description based on a Langevin-type equation could be introduced 
to gain information on fully quantum properties of simple linear open systems. In a forthcoming paper [3a | it will 
be shown that, by carefully dealing with the gauge freedom and the consequent dynamical constraints, the previous 
result can be extended to the case of N free quantum matter fields weakly interacting with the metric perturbations 
around a given background (here weakly interacting means that the gravitational coupling constant times the number 
of fields remains constant in the limit of large N). In particular, the correlation functions for the metric perturbations 
obtained using the Einstein-Langevin equation are equivalent to the leading order contribution in the large N limit 
to the correlation functions that would follow from a purely quantum field theory calculation. This will generalize 
the results already obtained on a Minkowski background [3a, |37[ . These results have important implications on the 
use of the Einstein-Langevin equation to address situations in which the background configuration for the scalar field 
vanishes, so that linearization around such a configuration is no longer possible. This includes not only the case of 
a Minkowski background spacetime, but also the remarkably interesting case of inflationary models driven by the 
vacuum polarization of a large number of conformal fields with vanishing expectation value [3ll . WA |88| , where the 
usual approaches based on the linearization of both the metric perturbations and the scalar field perturbations and 
their subsequent quantization can no longer be applied. 
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APPENDIX A: DEFINITION AND BASIC PROPERTIES OF GAUSSIAN STATES 

In this appendix we summarize the definition and basic properties of Gaussian pure states. In contrast to the rest 
of the paper, W will be used throughout this appendix to denote either the state of a field or its wave functional in 
the Schrodinger picture rather than a gauge invariant variable for scalar- type metric perturbations. 
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A pure state is called Gaussian if its wave functional in the Schrodinger picture is a Gaussian functional: 



$ [0(f)] oc exp 



d 3 xd 3 x'(j)(x)A(x,x')(t)(x') + / d 3 xB{x)<P(x) 



(Al) 



with a suitable normalization constant. The fundamental property of Gaussian states is the fact that the cumulants 
of order higher than two associated with the quantum expectation values of products of the field operator vanish, 



Sj(xi 



Sj(x n ) \ ' 



0, 



(A2) 



3=0 



for n > 3 and where we introduced the notation (O)* = (^lOI^). This implies that the connected part of any quantum 
correlation function 0(f„) . . . 0(fi) Y$>) with n > 3 vanishes or, equivalently, that any quantum correlation function 
(^| 0(f n ) ■ ■ ■ 0(f i) \*&) can be written as a linear combination of products involving the expectation values (^1 0(fi) \^} 
and two-point functions (^1 0(fj)0(ft) \^>). Furthermore, if the Hamiltonian of the field under consideration is 
quadratic, this property can be generalized for different times to quantum correlation functions in the Heisenberg 
picture of the form ^(^1 0(i„, f„) . . . 0(ti, x\) \^) H , which follows from Wick's theorem [891 ]. 

Finally, given a field operator and a Gaussian state jvP) with non- vanishing expectation value (</>)*, it is always 
possible to introduce a new field <f> = 0— (0)* so that the wave functional ^[(fi\ for the state \^>) in the basis associated 
with the field ip becomes a Gaussian functional with vanishing expectation value and independent of the expectation 
value (0)*. This can be immediately seen by rewriting the expression for the wave functional in Eq. (|A1[) as 



^ [0(f)] oc exp 



d 3 xd 3 . 



o' (0(f) - (0(f))*) A(x,x') (0(f) - (0(f))* 



and then change to the basis associated with the field (p: 



^> [(p(x )] oc exp 



d 3 xd 3 x'(p(x )A(x , x')ip{x') 



(A3) 



(A4) 



It is precisely in this sense that the state for the inflaton field perturbations (p introduced in Sec. Ill Bl follows imme- 
diately from the state of the inflaton field 0. 

It should be emphasized that any of the vacuum states commonly considered for free fields in curved spacetimes 
are Gaussian states. Furthermore, as stated in Sec. IIIB1 in this paper we concentrate on states which are invariant 
under the symmetries of the Robertson- Walker metric. In particular, for the case of Robertson- Walker metrics with 
flat spatial sections this implies (0)*(t, f) = (0)*(£) and A(x,x') — A(\x — f'|). 



APPENDIX B: CONVENTIONS AND NOTATION FOR THE LINEARIZED STRESS TENSOR 

Here we explain the notation concerning the linearized stress tensor employed in this article and clarify some related 
subtle points. 

Let us begin with the objects which appear in the Einstein-Langevin equation ([7]). The tensors and T^' 
correspond to the terms proportional to the metric perturbations in the perturbed version of the background (un- 
perturbed) objects G ®) and T^?'. The indices of G^) can be raised using the background metric [from now on 
everything that will be said for G^ applies exactly in the same way to T^}. On the other hand, one could perturb 
the background object G^ ab with the indices already raised, and reach a different result for G^ ah : it would differ 
by the terms G^ cd h^ + G^ ac h b c . Our notation (and those commonly employed) is ambiguous in the sense that it 
does not distinguish between both possibilities. In order to remove such an ambiguity it is necessary (and sufficient) 
to specify a priori which objects are going to be perturbed. In particular, in Sec. [TT] it is G^) and TjfP that are 
perturbed (of course, one can then freely use the background metric to raise and lower indices), whereas it is G a 
and T (0) a b that will be perturbed in Sees. |TTT] and Ml and finally G^ ab and T^ ab in Appendix O Fortunately, this 
will not change the form of the Einstein-Langevin equation since the terms corresponding to the difference between 
both choices for the Einstein tensor and the stress tensor cancel out because their background counterparts satisfy 
the Einstein equation = . Hence, if one deals with equations involving tensorial objects rather than with 
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isolated objects, everything is independent of the particular choice provided that the same choice is made for all the 
objects in the equation. 

The previous ambiguity does not affect the stochastic source of the Einstein-Langevin equation, which is completely 
defined on the background spacetime (the noise kernel is evaluated on the background geometry). Furthermore, the 
argument given in Sec. IIHI to show that $ = $ is not affected by such an ambiguity either. The reason is that the 
ambiguous extra terms for 5%j with i ^ j vanish because both (T$?) re n (or T^v) and the scalar metric perturbations 
in the longitudinal gauge are diagonal. Similarly, each one of the terms appearing in the Oi component of the Einstein- 
Langevin and the quantum version of the linearized Einstein equations considered in detail in Sees. IIII Bl and IIII CI 
respectively, do not suffer from the ambiguity either. 

The tensors an d , which are discussed above, result from linearizing just the metric perturbations. On 
the other hand, in the usual treatment of cosmological perturbations in inflationary models not only the metric 
perturbations, but also the inflaton perturbations are simultaneously linearized. Therefore, we introduced the notation 
T a b for the contribution to the expectation value of the background stress tensor (T^) Ten which is quadratic in the 
background solution 4>(t]) and independent of the inflaton perturbations, i.e., the first term on the right-hand side of 
Eq. (p~2|) : T ao = (T^)^. As we already explained in Sec. Ill Bl the expectation value {T^') Ien has also a contribution 
which is quadratic in the inflaton perturbations, but it is neglected when linearizing in those. We also introduced 
the notation 8T a b for the perturbed stress tensor operator obtained when linearizing with respect to both the metric 
perturbations and the inflaton perturbations. Its expectation value (ST a b[g + h]} is, thus, equivalent to linearizing also 
with respect to the inflaton perturbations the expectation value (T^l [g + h]) Tcn . Similarly, the operator 5T a b[g + h] 

corresponds to linearizing in the inflaton perturbations the expression Tjv — (Tf?) + . This expression, which 
appears (after raising one index) in Eq. (|26|) . might seem a bit awkward at that point, but this is just because the 
notation generally employed for the Einstein-Langevin equation, where one linearizes only with respect to the metric 
perturbations, is no longer the most natural when one also linearizes with respect to the inflaton perturbations. The 
expression is appropriate either when quantizing both the metric perturbations and the scalar field perturbations 
or when considering a stochastic version of it, namely the Einstein-Langevin equation. In the latter case one takes 
the expectation value of plus a stochastic source that accounts for the quantum fluctuations of the operator 
i a b = Tj$ — (T"*?), whose expectation value vanishes. Finally, the notation 6t a b is used for the result of linearizing 
the operator t a b with respect to the inflaton perturbations, which coincides with 8T a b evaluated on the unperturbed 
metric. 



APPENDIX C: FUNCTIONAL APPROACH TO THE EINSTEIN-LANGEVIN EQUATION AND 

ALTERNATIVE DERIVATION OF EQ. (f52|) 

The Einstein-Langevin equation for metric perturbations around agiven background and interacting with quantum 
matter fields has been formally derived using functional methods [53|, [54|, [55|, [56|, 157], [58|, l6fj • This was achieved by 
regarding the metric perturbations as an open quantum system with the environment corresponding to the quantum 
matter fields, and using the influence functional formalism for open quantum systems introduced by Feynman and 
Vernon (321 . |33| . In this appendix we briefly review some basic aspects of the functional approach to the Einstein- 
Langevin equation and explain how an alternative derivation of Eq. (|52p for the expectation value of the linearized 
stress tensor operator evaluated on the perturbed metric can be obtained. 

When considering derivations of the Einstein-Langevin equation using functional methods, one begins by computing 
the influence functional for the metric perturbations by integrating out the quantum matter fields (we will only consider 
free fields) as follows: 

rf{tf)=ipt pip' (t t )=<pi 

e iS lF [h,h\ = / dlpfd(p . dlp >. / Vlp V(p ' e iS[ V ,g+h]-iS[ V <9+h] p ^ (phip > i . t .^ (C1) 

J J <p(ti)=<pi Jip'(ti) — (p^ 

where p [tpi, tp[; t\) is the density matrix for the initial state of the matter field, tp, which is assumed to be initially 
uncorrelated with the metric perturbations (moreover, asymptotic initial conditions with ti — > — oo are usually con- 
sidered) and S [ip, g + h] is the action for the matter field evolving on a spacetime with metric g a b + h a b- Furthermore, 
only terms up to quadratic order in the metric perturbations h a b around the fixed background metric g a b will be 
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considered. In that case the action for the matter field can be written as 

1 



S[<p,g + h] = d^xyf-g (x)h ab (x) T ab [<p, g ab ; x) 

+ i / dixy/ ~ 9 ^ J d4 yV-9{y) h ab {%) Kd iy) 

1 8{^gJxjT ab [ V ,g ab -x)) 
y/-g{x)y/-g{y) Sg cd (y) 



O {h 3 ab ) , (C2) 



where T ab [tp, g' ab ; x) = 2(—g(x)) SS [ip, g] /8g ab (x) corresponds to the stress tensor for the matter field, whose 
functional derivative {—g(y)) 5T ab [(p,g ab ;x) /Sg c d(y) is a local object, i.e., proportional to the covariant delta 
function (— g (y)) 6^ (x — y). The influence action to quadratic order in the metric perturbations, which can be 
obtained by integrating out the matter field ip as explained in Refs. (53l. [60|. exhibits a structure analogous to that of 
a linear open quantum system: 

Sip [E, A] = Z ■ A + A • (H + K) ■ E + • N ■ A, (C3) 



where A ■ B denotes J d 4 x^—g (x)A ab (x) B ab (x) and we introduced the average and difference variables E (, = 
{h ab + h' ab )/2 and A a b — h' ab — h ab . The expressions for the kernels are the following: 

Z ab (x) = - l -(f ab [p l9ab -x)) : (C4) 

H abcd (x,y) - hm(T*f ab [p,g;x)f cd [p,g;y)) - l -([f ab [p,g;x),f cd [p,g;y)}), (C5) 



K^y) = - ^ ^a^ , (C6 ) 

2y/-g(x)y/-g(y) Sg cd (y) 

N abcd (x,y) = i<{£ afc [^ ff ;x),t crf [^, 5 ;y)}), (C7) 



where the functional derivative in Eq. (|C6[) should be understood to account only for the explicit dependence on the 
metric, whereas the implicit dependence through the field operator ip[g] is not included. The notation T* appearing 
in Eq. (|C5[) means that the matter field operators must be temporally ordered before applying any derivatives acting 
on them. Thus, we have, for instance, (T*V^(x)V^(y)) = Vj^V^ (T<p(x)(p(y)). Note that although the background 
geometry is non-trivial in general, the notion of temporal ordering is well defined because we are restricting the possible 
background geometries to globally hyperbolic spacetimes, which are time orientable; moreover, the microcausality 
condition of the quantum field theory for the matter fields under consideration guarantees that {[0\{x), 0%(y)]) = 
if Oi(x) and C>2(y) are local operators and x and y are spacelike separated points. It should also be noted that the 
first term on the right-hand side of Eq. (|C5p is symmetric under interchange of x and y, whereas the second one is 
completely antisymmetric. On the other hand, the term on the right-hand side of Eq. (|C6|) is local and symmetric 
under interchange of x and y. 

The noise kernel N abcd (x,y) requires no renormalization, as explained in Sec. Ill Al whereas the kernels Z ab (x), 
H abcd {x, y) and K abcd (x, y) contain diver gences [some regularization procedure is implicitly understood in Eqs. (|C4[) - 
(|C7[) ] that can be canceled out by adding suitable counterterms, quadratic in the curvature, to the bare gravitational 
action. These are precisely the same counterterms which arc introduced in scmiclassical gravity so that, when 
functionally differentiating with respect to the metric, they cancel the divergences from the expectation value of the 
stress tensor. This fact should not be surprising at all since the kernel Z ab {x) corresponds to the expectation value 
of the stress tensor operator on the background metric and the kernels H abcd (x, y) and K abcd (x, y) are closely related 
to the expectation value of the stress tensor operator on the perturbed metric, as follows straightforwardly from the 
following relation, valid up to linear order in h ab : 



(f" b [g + h;x))= 2 *S*[h,h1 



v / -(g + /i)(x) Shab(x) 



-2Z ab {x) - 2(H ■ h) ab (x) - 2(M ■ h) ab (x), (C8) 



h'=h 



where we introduced the kernel M abcd (x,y) defined as follows: 

-1 S 
V-5(y) s 9cd(v) 



M^(x, y) , - 1 (C9 ) 
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which results from adding to the kernel K abcd (x,y) the term — (— g(x)) 1 ^ 2 {5\/~g(x)/8g c d{y))Z ab coming from the 
contribution to the factor 2[— (g + /i)(cc)] -1 / 2 that is linear in h ab . When the counterterms introduced in the bare 
gravitational action are included in the influence action, so that the divergences cancel out and the bare kernels 
H abcd (x,y) and M abcd {x,y) get renormalized, Eq. (|C8|) becomes 

= (f ab [g; x)) lcn - 2 (H Icn ■ h) ab (x) - 2 (M ron • h) ah (x), (CIO) 

h'=h 

which can be rewritten as 

(f^ ab [g + h; x)) Icn = -2 (H Icn ■ h) ab (x) - 2 (M ron • h) ah (x), (Cll) 

where, as mentioned above, the functional derivative appearing in the kernel M abcd (x,y) should be understood to 
account only for the explicit dependence of the stress tensor on the metric, whereas the implicit dependence through 
the field operator ip[g] is entirely contained in the first term on the right-hand side of Eq. (|C11[) . 

Taking into account the previous results, the Einstein-Langevin equation can then be obtained from the CTP 
effective action for the metric perturbations by using a formal trick described below. Such a CTP effective action 
for the metric perturbations has the following form at tree level (note, however, that the matter fields, which have 
already been integrated out, were treated beyond the tree level): 

r<?> p [h, h'\ = s e [h] - s g [h>] + 4 cn) [h, h>] + o (h 3 ab ) , (ci2) 

where S g [h] is the Einstein-Hilbert action J d 4 Xy/^gR up to quadratic order in the metric perturbations, and the finite 

parts of the local counterterms have been included in Sj^ n K On the other hand, using the following mathematical 
identity for the imaginary part of the influence action: 

e -im&p = e -iA-iv.A = det ( 27 rN)-^ /^e-^ ^ 1 ^-^, (C13) 

and interpreting £ ab as a stochastic source with vanishing expectation value and correlation function (t; ab (x)£ cd (y)) ^ = 
N abcd (x,y), one can define a stochastic effective action, 

r stoch [E, A] = S s [h] - S e [h'\ + A • (H rcn + M rcn ) . E - If . A, (C14) 

such that (exp(ir st0 ch))^ = exp(ir[^! r p). The Einstein-Langevin equation can be immediately obtained by functionally 
differentiating with respect to the metric perturbation h a b and letting h' ab — h ab afterwards: 

= _I G (D ^ [g + h . x) _ . h) ob {x) _ (Mrcn . h) ab {x) + Kab {x) ^ 

Zk 2 

h'=h 

It is worth discussing the issue of causality in Eq. (|C15[) . which basically amounts to considering the second term 
on the right-hand side of the equation since the remaining terms are local. The right-hand side of Eq. (|C5|) can be 
formally rewritten as 

- ~ ( [t ab [0, g; x) , t cd [0, g; y)} ) 9* ( Vx - r, y ) , (C16) 

where r\ x and r\ y can be any pair of well-behaved time coordinates for the points x and y, and the star index in the 
theta function was used to indicate that the derivative operators acting on the scalar field which appear in the stress 
tensor operator should also act on the theta function. Thus, all the terms in expression (|C16[) are either proportional 
to 0(r] x — r] y ), S(r] x — rj y ) or S'(j] x — r/y), and, being proportional to a commutator, expression (|C16[) vanishes for 
spacelike separated points because of the microcausality condition of the quantum field theory for the matter fields. 
Furthermore, since both the divergences and the counterterms are local (proportional to delta functions or derivatives 
of them), the contribution to Eq. (|C15p from the term (H Yen ■ h) ab (x) is causal, i.e., it only depends on the metric 
perturbations h c d(y) at any point within the past lightcone of x. 
Finally, taking into account Eq. (|C11[) . Eq. (|C15|) becomes 

G { 2 \9 + h] = K(f£> [g + h]) rcn + < ah , (C17) 



(f ab [g + h;x)) Icn = 



SS 



(rcn) 
IF 



[h, h'} 



^-{g + h)(x) Sh ab (x) 



q _ i <5r st0 ch 

y/-(g + h){x) Sh ab (x) 
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where the indices have been lowered using the background metric. It should be noted that, in contrast to Sec. Ill A[ 
the tensors appearing in Eq. (|C17|) correspond to perturb the background tensors with both indices already raised. 
However, as pointed out in Appendix [Bj the resulting equations in both cases are equivalent because the unperturbed 
tensors satisfy the semiclassical Einstein equation. Therefore, Eq. (|C17[) is in complete agreement with Eq. |(7J), 
keeping in mind that the finite contributions of the counterterms, corresponding to A a b and B a b in Eq. ([4]), have been 
reabsorbed in the renormalized expectation value of the stress tensor operator. 

After this brief review of the functional approach to the Einstein-Langevin equation, let us now see how Eq. (|C11|) 
gives a result for (67q)& which is equivalent to that obtained in Sec. lIIIBl To begin with, it should be pointed out that 
the ambiguity mentioned in Appendix [B] does not affect the Oi component of (ST^)® since both the background stress 
tensor and the scalar metric perturbations in the longitudinal gauge are diagonal. Furthermore, it can be seen that 
for a diagonal perturbed metric the Oi component of the second term on the right-hand side of Eq. (|C11[) vanishes. 
Thus, we can concentrate on the first term. 

Fourier transforming the spatial coordinates as done in Sec. IIII1 the expectation value for the Oi component of the 
perturbed stress tensor becomes 

((T^)k[g + h])(rj) = -2 1 ^/^W)dr)'H 0icd (v,ri'; k)h c k d (r,'), (C18) 

where the kernel H a f, c( i(T},r)';k) corresponds to the Fourier transform of the two terms on the right-hand side of 
Eq. (|C5|) . Using the equivalent expression in Eq. (|C16[) . the kernel H a bcd(i], T)'', k) is given by the following expression, 
which already takes into account that the Fourier transform of the expectation value / [ia&[<?; ??, x), i cd [g; if, ) is 
proportional to a Dirac delta function due to the existing translation invariance in the spatial coordinates, 

H abcd (v,ri';k)(2TT) 3 5(k + k') = -\([{i ab ) k [g;v), {ic d ) k , [g;v')]) 0*(v-v'), (C19) 

where the star in the theta function had been introduced earlier to indicate that the derivatives appearing in t a b and 
t c d should also act on the theta function. Performing a similar decomposition to that introduced for the noise kernel 
in Eq. (|20p . we obtain two non-vanishing contributions to the expectation value ( \t a b\x; x )> tcd[g', x')] ) Iff]- 

< [iab[g; x),icd[g; x')] ) = ( [i a b[g; x), t C d[g; x')] ) <j)tf> + ( [i ab [g; x)j cd [g; x')] ) w , (C20) 

where the first contribution is quadratic in the quantum operator <p[g] for the inflaton perturbations evolving on the 
unperturbed geometry, whereas the second contribution is quartic in <p[g]. As already pointed out for the separation 
of the noise kernel, the fact that the conservation of the stress tensor, which is the source of the Einstein equation, is 
necessary to guarantee its integrability implies that both contributions to the expectation value must be separately 
conserved if we want to discard one of them keeping the consistency of the Einstein equation at the order that we are 
working, which is linear in the metric perturbations. This is indeed the case as follows from the fact that both the 
background homogeneous solution 4>{rj) and the operator (p[g\ satisfy the Klein-Gordon equation on the background 
spacetime. 

If we keep only the first term on the right-hand side of Eq. (|C20|) . i.e., if we take H a b c d(x,x') = 
( \8tab[g') x), St c d[g; x')] ) 9*{r] x — r)' x ), which corresponds to considering the contributions to t a h[g] that are linear in 
the inflaton perturbations and is consistent with the linearization of the inflaton perturbations that was considered 
in Sees. IHII and IIV1 we obtain 

((Sf 0i ) k (v)h = -i<t>'mh) f dr,'a 2 ( V ')<i> k (v') {W) ([0k fa), <P'- k (v')}) 

J no 

-2mV(7 ? ')^(r?')([^fa)^-fcfa / )])}^fa^') ) (C21) 

where we used the explicit expressions for the components -f/oioo an< i Hoijj- I n this case, there is actually no need for 
the * prescription in 9*{r]^rj'), which implies that the derivative acting on <p{rf) should also act on the theta function, 
since it yields a term proportional to [ip(j]),ip(j])]d(r] — r]'), which vanishes identically. Integrating by parts the first 
term in the integrand and using the Klein-Gordon equation for the background solution 4>[rf), given by Eq. (|17p . we 
finally get 

((SfjUv))* = ia-*{r,)<t>'{rj)(ih) f d v >a 2 ( v >) {4[Mv), 0-k(v'W(v'Wk(v') 

J no 

-2[0 k {n), 0_ k {r ] ')]m 2 a 2 {r l ')4>(v')^kiri')} , (C22) 
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where the factor a~ 2 (rj) comes from raising the index i with the background metric and we have substituted the 
expectation value for the commutator of the field operators simply by the commutator since for a linear theory they 
are c-numbers, whose expectation value is independent of the state. This result for the expectation value of the stress 
tensor coincides with Eq. (|52"|) , found in Sec. IIIIBI It should be noted that the contribution from the boundary term 
at rj' = 77 which results from the integration by parts vanishes because [ip{r]), <p{r])] = 0. On the other hand, there is 
a non- vanishing contribution from the boundary term at rj = 770: 



It might seem that the existence of this term would imply a conflict between the result for the expectation value of the 
linearized stress tensor operator obtained in Sec. IIII 51 using the equations of motion for the quantum operators in the 
Heisenberg picture and the result based on the influence functional formalism derived in this appendix. However, this 
is not the case. The reason for the apparent discrepancy is the following. When computing the expectation value of the 
stress tensor operator, there are terms proportional to (<p\g; 7])}, where the operator <p\g; 77), which satisfies the Klein- 
Gordon equation, can be written as a linear combination of a term proportional to (p[g; 770) and a term proportional 
to (p'\g; r)o). In particular, in Sec. IIIII we chose a state for which both (<p[g',r]o)) and (<p'[g; 770)) vanished. On the 
other hand, in the approach based on the influence functional formalism the operators which naturally determine 
ip[g; 77) in terms of the initial state are ip[g; rjo) and its conjugate momentum n[g; 770). Since the coupling between the 
metric perturbations and the inflaton perturbations involves terms proportional to the time derivative of the inflaton 
perturbations, fc[g; 770) will differ from <p'[g; r/o) by a term proportional to the metric perturbations at the initial time. 
This is precisely the origin of the term in expression (|C23|) . Thus, the apparent discrepancy is just a consequence of 
the fact that in the influence functional approach it has been implicitly assumed that the initial state has vanishing 
(^[5; 770)) rather than vanishing (<p'[g; 770)) - 

It is important to stress that the expression in Eq. (|C22|) for the expectation value of the stress tensor operator needs 
no renormalization. This fact can be easily understood because we are dealing with the linearized theory. Therefore, 
the terms involved in the computation of the expectation value of the stress tensor operator are proportional to 
(ip[g + h; x]}, whereas the divergences that arise in an exact treatment (without linearizing with respect to the scalar 
field) are a consequence of taking the coincidence limit x — > x in terms involving products of the field operator, i.e., 
proportional to (ip[g + h; x]ip[g + h; x 1 ]). Alternatively, when considering Eq. (|C8|) together with Eq. (|C16[) . the need 
for renormalization can be understood as follows. The expectation value of the commutator is finite as long as one 
restricts to x ^ x', but it diverges when one considers the coincidence limit. Nevertheless, it is still meaningful as 
a distribution. In this context, the divergences arise because the product of distributions in Eq. (|CJ16[) is ill defined 
in general although each factor is well defined as a distribution; see Ref. [90( for a detailed discussion on this point. 
In fact, the terms in Eq. (|C16[) involve terms proportional to the imaginary part of the product of two Feynman 
propagators Gp(x, x')Gp(x, x') [13, HH, [2H . Working in Fourier space for the spatial variables, this product becomes 
J d 3 qG f (77 , rj '; k — c[)Gf (77, 77'; q), which exhibits an ultraviolet divergence when performing the integral J d 3 q over 
all possible momenta. On the other hand, when linearizing with respect to the scalar field, the Fourier transformed 
version of the terms in Eq. (|C16p is simply proportional to £^(77,77'; fc), with no integral over momenta and, hence, 
no ultraviolet divergence. 



In this appendix we will explain why, strictly speaking, a homogeneous solution & k (77) with some particular initial 

conditions must be added to the purely inhomogeneous solution (77) when solving the Einstein-Langevin equation 
as done in Sees. [TTT1 and HVl The situation is completely analogous when solving the quantum version of the linearized 
Einstein equation as in Sec. IIII CI 

It is well known that the Bianchi identity guarantees the integrability of the Einstein equation provided that the 
stress tensor of the matter sources is covariantly conserved. Let us, however, discuss this point in some more detail. 
The ten components of the Einstein equation for a globally hyperbolic spacetime, which can be foliated with a set 
of Cauchy hypersurfaces, can be formulated as an initial value problem with time corresponding to some continuous 
variable labeling the Cauchy hypersurfaces. In particular for the cosmological problem that we are considering we can 
choose the homogeneous spatial sections labeled by the conformal time 77 as the set of Cauchy hypersurfaces. The four 
temporal components of the Einstein equation can then be regarded as a set of dynamical constraints at any given 
instant of time. Thus, the integrability of the Einstein equation as an initial value problem can be understood in the 
following way: using the Bianchi identity and the conservation of the matter sources, the constraints can be shown 
to hold at any time provided that the spatial components of the Einstein equation are satisfied for all times and the 




(C23) 
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four constraint equations are fulfilled on the Cauchy hypersurface corresponding to some initial time [40( • Obviously, 
the previous discussion can be extended to the case of the Einstein-Langevin equation since the stochastic source is 
also covariantly conserved. 

Let us recall the temporal components of the Einstein-Langevin equation for scalar metric perturbations after 
Fourier transforming with respect to the spatial coordinates: 

(k 2 + m 2 ) <£ k + m& k = V (< W)*>* + (£ °) fe ) , (di) 

iki@' k + H$ k ) = |a 2 (((<5T l ) fc ) $ + (C ) k ) . (D2) 

In Sees. UlTl and HVl the constraint equation (|D2[) was solved to find ^ k (rj). However, one should make sure that 
the remaining components of the Einstein-Langevin equation are also satisfied. According to the discussion in the 
previous paragraph, to make sure that this is indeed the case it is sufficient to demand that the Eq. (|D1[) holds at 



the initial time 770 for every k. The solution of Eq. (|D2p can always be written as ^> k (r]) = <& k (rf) + & k (?]), where 

3>j[ (77) is a solution of the inhomogeneous equation which vanishes at the initial time 770 and $1 (77) is a solution of 
the homogenous equation which is completely determined by specifying its value at 770. Imposing Eq. (|D1|) at 770 and 

(h). 



using Eq. (|D2[) evaluated at 770 in order to substitute &' k {T]o), one obtains the following result for <j>^ (771 



(D3) 



where we took into account that ((57q )fc(?7o))# vanishes, as can be immediately seen from Eqs. (|52|) . or (|C21|) . because 
the limits of integration coincide. We also used the fact that ((#?o )fe(»7o))* = a ^ 2 ( r lo)W( r lo)] 2 ^k('no)' in this case 
the second term on the right-hand side of Eq. (|C8[) vanishes for the same reason as with {(STq )fc(?7o))$! but there is 
a non- vanishing contribution from the last term in Eq. (|C8[) . which corresponds to the first term on the right-hand 
side of Eq. ([29]). 

Since Eq. (|D2|) is a first order integro-differential equation, the result for & k h \r/o) in Eq. (|D3|) completely determines 

the homogeneous solution & k h ^ (77) . The situation will be completely analogous when linearizing and quantizing both 
the metric perturbations and the inflaton perturbations, as done in Sec. IIII C[ with the quantum operator for the 
metric perturbations $(2;) replacing the stochastic scalar field &(x) and the operator 5t h a instead of the stochastic 
source £„(x). Hence, the argument concerning the equivalence between the quantum correlation function for the 
metric perturbations and the stochastic correlation function can be straightforwardly extended, following the same 
line of reasoning as in Sec. IIII C\ to the case in which the contribution from the homogeneous solution is also taken 
into account. 

Nevertheless, in Sees. Mil and IIVI this homogeneous solution was not considered when giving the final result for 
the correlation function of the metric perturbations. Therefore, we end this appendix arguing why it is justified to 
neglect the contribution from the homogeneous solution when computing the correlation functions for scalar metric 
perturbations at large scales in the context of cosmological inflationary models. In other words, the contribution 
from the first three terms on the right-hand side of Eq. (|45p is much smaller than the contribution from the fourth 
term when considering a situation similar to that addressed in Sec. IIVI This can be qualitatively understood in the 
following way. Since Eq. (|D2|) is a linear first order differential equation, the solution of the homogeneous equation, 



$k (77), will be proportional to the expression for § { k \rj) given by Eq. (|D3|) . Thus, the first term on the right-hand 
side of Eq. (|45f is proportional to the correlation functions for the stochastic source at the time 770, and the second 
and third terms are proportional to the correlation functions at different times: 770 and a time 77' which is integrated 
from 770 to 771 or 772 [see Eq. f|44[) ] . The value of the noise kernel is small when one or both the two arguments are 770 
provided that 770 is negative enough so that the scales of interest were well inside the horizon at that time. This is in 
contrast to the contribution from the last term in Eq. ([45|) when the relevant scales are well outside the horizon at 771 
and 772, since the two arguments of the noise kernel in that term are integrated from 770 to 771 or 772. Hence, the reason 
for neglecting the first three terms on the right-hand side of Eq. ([IS")) in this context is actually rather similar to the 
reason for the weak dependence on 770 of the result obtained in Sec. IIVI when the scale k is well outside the horizon at 
771 and 772, and 770 is negative enough so that k is well inside the horizon at that time. 

The previous argument can be made more precise if we concentrate on the particular model considered in Sec. IIVI 
In that case, if we neglect the non-local term corresponding to (57q)$, as done in Sec. IIVI the expression for the 
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homogenous solution is 



2 a( V ) 



3>H(Vo 

i k; 



(D4) 



The contribution form the first three terms on the right-hand side of Eq. (|45l) can then be explicitly computed and 
compared to the last term, taking into account that krji, kr\\ <C 1 and ki]Q 1. In particular, the first term on the 
right-hand side of Eq. (|45|) is proportional to (K/2)(m/m p ) 2 (2Tr) 3 k~ 3 S(k — k')a 2 (r]o) / a(T]i)a(r]2) and a sum of terms 
of order 1, (l/fc?7o) 2 and (m/H)(l/kr]o) 2 . The factor a 2 (r]o)/a(r]i)a(r]2), which is of order krnkr]2 / (krjo) 2 , as well as 
(1/krjo) 2 and m/ff are much smaller than 1. It is thus clear that those contributions can be safely neglected as 
compared to the last term, which was found to be of order (K/2)(m/m p ) 2 (2ir) 3 k~ 3 5(k — k') in Sec. IIV1 

Similarly, the second and third terms on the right-hand side of Eq. (|43)l are proportional to 
(K/2)(m/m p ) 2 (2,Tr) 3 k~ 3 6(k — k')a 2 (rj ) / a(rji)a(i]2) and a sum of terms of order 1 and 1/krjo. Therefore, they can 
also be neglected as compared to the last term. 



APPENDIX E: ALTERNATIVE PROOF OF THE EQUIVALENCE BETWEEN STOCHASTIC AND 

QUANTUM CORRELATION FUNCTIONS 

In this appendix we provide an alternative proof of the equivalence between stochastic and quantum correlation 
functions whose key step is to show that the Einstein-Langevin equation for linearized cosmological perturbations 
implies Eq. (6.48) of Ref. [z3- 

Let us consider the Einstein-Langevin equation for scalar metric perturbations when one also linearizes with respect 
to the inflaton field, whose different components are given by Eqs. (|5T)) - ([55)) . We will take '5 = <&, as justified by the 
discussion before Eq. (|40p . and work in Fourier space for the spatial coordinates. Next, we add Eq. (|37|) . the i = j 
component of Eq. and Eq. multiplied by F = 2m 2 a 2 ((f>/(f>')(ki/k 2 ), which leads to the following result: 



I + 2 (h - *' fc + k 2 <P k +2 Oh'- H^j $ fc = \a 2 [(5%°)* + (ffZ?)* + F(5T % + £ ° + & + F%] ^ , (El) 
with no summation over the repeated i indices. In deriving Eq. (|E 1 1) we made use of the following two relations 

6W + 2m 2 a 2 -^ = % (m(/)' + m 2 a 2 <j)) = 2H-2^-, (E2) 

q/ qy q> 

2H' + 2H [2H + m 2 a 2 ^- ) =2H'-2H^- r , (E3) 
V 9 ) <? 

which follow from the Klein-Gordon equation (jTTJ) for the background field <j>. The final step is to show that the 
right-hand side of Eq. (|E1[) vanishes. In order to do so, it is convenient to consider first the Fourier-transformed 
version of Eqs. (f2l)|) - PTj) for the linearized stress tensor. It is then straightforward to show that 

(S^) k + (SV) k + F(S^) k = 0, (E4) 

with no summation over the repeated i indices. The same conclusion applies when tp is promoted to a Heisenberg oper- 
ator ip, which implies that the first three terms on the right-hand side of Eq. (|E 1 1) cancel out. On the other hand, since 
(£o + £| + F£,q) k is a Gaussian stochastic process with vanishing mean, in order to prove that it vanishes it is sufficient 

to see that (($ + Q + F$) k (ri) (^) fc , (rf)) , which is proportional to ({(«Jt° + St\ + F5i ) k ( V ), (5%) k , (r/')}), van- 

ishes. Indeed, taking $ = in Eq. (|E4|) and promoting (p to a Heisenberg operator, it follows that (5t$ + 5t\ + FSt l ) k = 
0. Thus, the right-hand side of Eq. (|E1[) vanishes and one is left with 



+ 2\H- ^-J - V 2 $ + 2 \H' - n ^Fj ® = (E5) 

which coincides with Eq. (6.48) in Ref. [74|. Several remarks about Eq. (|E5|) are in order. First, the non-local terms 
associated with (ST^}® are not present so that, when working in Fourier space for the spatial coordinates, one is 
left with an ordinary differential equation rather than an integro-differetnial one. Second, the equation exhibits no 
dependence on the stochastic source. However, the solutions of the Einstein-Langevin equation should also satisfy 
the constraint equations at the initial time in addition to Eq. (|E5[) . According to the results in Appendix [D] this 
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implies a dependence on the stochastic source for the initial conditions $k(vo) an d ^'kivo): which will involve a linear 
combination of terms linearly proportional to the stochastic source [as given by Eq. (|D3p and an analogous result for 
$' fc (r/ ) that can be obtained by substituting Eq. (|D3j) into Eq. (|D2j) ]. The solution of the linearized Einstein-Langevin 
equation can then be written as Qk(v) = u i(v)^k(^(vo)]Vo)+ u 2(v)^'k(£.'il( r lo)] r lo)i where ui{rf) and 112(77) are solutions 
of Eq. (|E5[) with initial conditions Mi(^o) = 1> ^i(?7o) — and 1*2(^0) = 0, u' 2 {r]o) = 1. Such a dependence of the 
initial conditions on the stochastic source at the initial time is responsible for the non-trivial part of the stochastic 
correlation functions of <& at later times. 

If one quantizes both the linearized metric perturbation and the inflaton field, so that $ and ip are promoted to 
Heisenberg operators in Eqs. ([33|) - (|35)) . one can easily conclude (proceeding analogously to the previous paragraph) 
that the operator $ for the linearized metric perturbations also satisfies Eq. (|E5jl . Furthermore, the constraints at 
the initial time give the same results for &k(vo) and ^(rjo) as in the stochastic case but with St" in place of 
so that the solutions of the quantum version of Eqs. (f33|) - (j35|) can be written as 'tkiv) — u i( J 7)^fc(<^( J 7o); Vo) + 
MvWk(K( m ),m)- Taking into account that (fy [<?; x) £ [g; y)) $ = (1/2) ({6t» [g; x) , 6% [<?; »)}), it is straight- 

forward to see that the result for symmetritzed two-point quantum correlation function (l/2)({$fc(r7i), $£'(772)}) is 
equivalent to that for the stochastic correlation function ($k(vi)$k' (%)){ • This constitutes an alternative proof to 
that provided in Sec. IIII CI of the equivalence between quantum and stochastic correlation functions for linearized 
cosmological perturbations. 

We close this appendix by discussing a recent claim that there is a discrepancy between stochastic gravity and 
the standard treatment for superhorizon modes. More specifically, in Ref. [93j the Einstein-Langevin equation for 
linearized cosmological perturbations was solved using certain approximations and the correlation function for the 
gauge invariant variable C was computed (this variable corresponds to the curvature perturbation in the uniform 
density gauge 94], which coincides with the comoving gauge for modes outside the horizon [§5j|). It was found that 
for modes outside the horizon and for a sufficiently large number of e-folds the correlation function was not constant 
in time, contrary to the standard result. This would be in conflict with the equivalence for linear perturbations that 
we have shown to hold in general. However, one can provide an exact argument which shows that the linearized 
Einstein-Langevin equation actually implies that for modes outside the horizon the correlation function of £ remains 
constant in time. This means that the result in Ref. [93[ seems to imply a problem with some of their approximations 
rather than a shortcoming of stochastic gravity for that regime as concluded there. 

The exact argument is the following. One starts with the expression for £ in terms of $ (z!, [§S H^l : 

C = 5 ^(* + W -V) + * I (E6) 

with w = p/ p, where the relation between the background density p and pressure p and the background field <j) is 
given by Eqs. (|13[) - (fLT| and the text after them. Differentiating Eq. (|E6[) with respect to the conformal time and 
multiplying by (3/2)W(l + w), one obtains 

3(1 l W) H£ = $" + 2 - ^ W + 2 (W - H^pj *, (E7) 

where we made use of the following relations, which can be derived from Eqs. (|15[) -([T 



l + w - — -^pr- — ?7j2 — . y^ 8 ) 



p + P = njcj)') 2 = 2(H 2 - H') 

p m 2 m 2 

£ Hl+v , )= 2$-2%. (E9) 



Note that the right hand-side of Eq. (|E7[) coincides with the left-hand side of Eq. I|E5[) except for the V 2 $ term. This 
means that for modes outside the horizon, for which V 2 $ can be neglected, the right-hand side of Eq. (|E7[) vanishes 
and the mode remains constant in time. 
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